NASA Contractor Report 3692 


Dynamic Effects of Internal 
Spur Gear Drives 


Adam Pintz, R. Kasuba, 

J. L. Prater, and R. August 


GRANT NAG3-186 
JUNE 1983 




NASA Contractor Report 3692 


Dynamic Effects of Internal 
Spur Gear Drives 


Adam Pintz, R. Kasuba, 

J. L. Prater, and R. August 

Cleveland State University 
Cleveland, Ohio 


Prepared for 

Lewis Research Center 

under Grant NAG3*186 


ms!\ 

National Aeronautics 
and Space Administration 

Scientific and Technicai 
Information Branch 



TABLE OF CONTENTS 


CHAPTER page 

I. INTRODUCTION 1 

1.1 General Remarks 1 

II. LITERATURE REVIEW 9 

2.1 External Spur Gears 9 

2.2 Internal Spur Gears 25 

III. ANALYTICAL INVESTIGATION 30 

3.1 Problem Pormulation 30 

3.2 Assximptions 35 

3.3 Method of Solution 37 

3.4 Static Analysis 38 

3.4.1 Nomenclature 38 

3.4.2 Local and Global Coordinate 

Systems 39 

3.4.3 External and Internal Gear 

Tooth Profile 41 

3.4.4 Contact Points Between Gear 

Tooth Pairs 48 

3.4.5 Line of Action, Contact Ratio 

and Interference Conditions 54 

3.4.6 Deflections and Stiffiiess of 

the Teeth and Their Supporting 
Structure 56 

3.5 Dynamic Analysis 62 

3.6 Computer Program 74 

3.6.1 Program Structure 74 

3.6.2 Description of the Executive 

Programs and Subroutines 77 

3.6. 2.1 Module 1 ^7 

iii 



CHAPTER 


PAGE 


3-6,2.2 Module 2 78 

3.6.2. 3 Module 3 79 

IV. RESULTS, DISCUSSION AND SUMMARY 80 

4.1 Results and Discussion 80 

4.1.1 Introduction 80 

4.1.2 Static Analysis 81 

4. 1.2.1 Comparison of ISG and ESG Set 

Perfonticince 81 

4. 1.2. 2 ISG Drives of Practiced. Interest . . 88 

4.1.2. 3 Effects of Radial Deflection on 

Static Performance 89 

4.1.3 Dynamic Analysis 91 

4.1. 3.1 Comparison of ISG and ESG Drive 

Dynamic Performance 91 

4.1. 3.2 ISG Drives of Practical Interest . . 102 

4.1. 3.3 Effect of Radial Deflection on 

Dynaunic Performance 102 

4.2 Summary and Conclusions 105 

BIBLIOGRAPHY 108 

APPENDICES 110 

A Spur Gear Formulae cind Involute Profile 

Development 110 

A-1 Standard Spur Gear Relations for 

the ISG Drive Ill 

A-2 Development of the Involute Profile . . 114 

B Deflections 118 

B-1 Bearing Deflections 119 

B-2 Radial Ring Gear Deflection 120 

B-3 Circumferential Deformation of 

Gear Teeth 121 

iv 


PAGE 


APPENDICES (Cont’d) 

B-3.1 Deflection of Point of Contact 

Due to Deformation of Teeth 121 

B-3,2 Deformation of the Teeth Due to 

Rotation of Their Foundation 125 

B-3.3 Deflection of the Teeth Due to 
Circumferential Deformation of 
the Rim and Gear Ring 128 

B-3.4 Hertzi2m Deformation at Contact 

Point 1 32 

C Computer Program Package 135 

C-1 Listing and Sample Run of the 

Static Analysis Program "Internal 

Static" 136 

C-2 Listing and Sample Run of the 

Dynamic Analysis Program "Internal 
Dynamic" 198 

C-3 Listing and Sample Run of the Stress 

Analysis "Intemad Stress" 248 

C-4 Entering of Input Data 271 

D Glossary of Terms 276 


V 


1 


CHAPTER I 
INTRODUCTION 


1,1 GENERAL REMARKS 

Sp\ar gears have been utilized for many years and are of great 
in^rtance in trauismitting power from one rotating shaft to another* 

To affect power transmission, two or more gears are combined in a 
variety of arrangements. The most conation and best understood config- 
uration is the use of two external spur gears side-by-side as shown 
in Figure 1. This arrangement of the gears is used in single and 
miiLtiple pairs or stages, and is referred to as a pcurallel shaft gear 
transmission or **gear box”. Parallel shaft gear boxes are economic 
in the power rzmge 0 to 1500 KW, but become large, heavy and less 
economic above these power levels. 

More compact arrangements are achieved with the use of an external 
gear inside an internal gear as shown in Figure 2. This configuration 
is referred to as internal gear drive cuid is applied for the movement 
of turntables, tank turrets, radar systems, and the transmission of 
power in wind turbines, wheel drives of off-the-road vehicles, etc. 

The internal gear drives can be used either as speed reducers or speed 
increasers . 

Because of their inherently higher cost, selection of internal 
gear drives is predicated on the need for compactness or the opportunity 
of sharing components with another f\anction in the system. Examples 
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of either possibility are shown in Figures 3 and 4. rhe radar gear 
reduction unit of Figure 3 is actually a combination or hybrid of 
external and internal gear reduction stages. The compactness of the 
design is obvious. 

Figure 4 is an example where two functions are combined in one 
component. In this application the rotor of a wind turbine is sup- 
ported at the inner race of a large ball bearing. Also, the inner 
race is used as the first stage speed increaser for the wind turbine. 

As shown in Figure 4, internal gear teeth at the bore of the inner 
race engage with the external mating gear. The mating gear is mounted 
at the far side of the beeuring. A long drive shaft connects the first 
stage speed increaser to the two-stage final speed increaser. Finally, 
the electric generator is connected to the final speed increaser and 
is shown at the far side of the drive shaft. The cost of a separate 
bearing support for the rotor and a three— stage speed increaser is 
higher than the arrangement shown in Figure 4. 

Other examples of internal gear arrangements can be found in 
epicyclic gearing (Figure 5) . As can be seen, an epicyclic gear train 
has a central "sun" gear, several "planets" meshing with the sun and 
spaced uniformly aroiind the sun, and an internal gear or ring gear 
meshing with the planets. The name epicyclic is derived from the fact 
that points on planets trace out epicycloidal curves in space . Because 
of the multiple use of plemet gears, epicyclic gearing is the most 
compact arrangement of all the spur gear systems. 

At present, concerted efforts are being made to increase the 
power-to-transmission weight ratio. These efforts are not limited to 
the aerospace industry alone because of the high cost of material. 
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Figure 3 - Assembly Drawing of Radar Gear - Reduction Unit 
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labor and energy. Savings in weight normally imply less material with 
some savings in cost of the transmission. 

Consequently, extension of the state-of-the-art in spur gearing 
is a continuing requirement. Unfortunately, spur gearing investiga- 
tions cure concentrated almost exclusively to external gears as is 
evident from the literature review in the next section. For example, 
the latest investigations utilize large scale digital computer pro- 
grams to analyze external spur gearing for both static and dynamic 
conditions. These sophisticated computer programs for external geaurs 
help to move the analytical simulation closer to the actual behavior. 
Prior to this research work, no such analytical tools were known to 
exist for internal spur gear (ISG) drives. 

Bie design of ISG drives was based principally on extrapolations 
of external gearing procedures by the American Gear Manufacturer's 
Association (AQIA) and the International Standard Organization (ISO) . 
In either case, these procedures represent the technology known in 
the 1950' s. The basic drawback of these procedures is that they are 
beised on highly idealized relationships, which in real applications 
hardly exist. TSius, a multitude of saifety and application factors 
are imposed on the procedures which can result in considerable 
overdesign. 



CHAPTER II 


L ITERATURE REVIEW 


As mentioned earlier the study of spur gearing is concentrated 
to a great extent on investigations of external spur gears. There- 
fore, this literature review relies heavily on the information for 
external spur gears in showing the progress and current status of 
spur gear technology. The information is presented chronologic cdly, 
and in separate sections for the external and internal spur gearing 
respectively. 

2.1 EXTERNAL SPUR GEARS 

In 1892 Lewis used the form of the gear tooth as one of the 
factors in a formula for the Strength of External Gear Teeth . He 
related the tooth load to the material working stress using simple 
beam theory cmd developed equation (1) , known as the Lewis Equation 
W « SPFY ... (1) 

vrtiere 

W = transmitted load, lb. 

S = safe working stress in material, psi 
P = circular pitch, in. 

F = face width of gears, in. 

Y * tooth- form factor 


♦Superscripts refer to entries in references. 
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Lewis also recognized that the instantaneous load of the teeth was 
affected by the velocity of the system. Barth took note of this 
fact and developed a formula which resulted in an adjustment of the 
allowable stress as follows: 


S 


d 


600 

600 + V 


...( 2 ) 


v^ere 

S ■ safe static stress, psi 
V * pitch line velocity, fpm 

This modified design stress was then used as the design stress 
in the Lewis Equation. Today, the American Gear Manufacturer’s Associ- 
ation (AGMA) reconmended practice for bending strength uses the Lewis 
Equation in modified form. 

In the 1920 's and the early 1930's the Americeui Society of 
Mechanical Engineers (ASME) Research Committee investigated gear tooth 
loads and available design criteria in order to develop a unified 
approach to gear design. Tests were conducted by Lewis auid Buckingham 
to determine the effects of production errors and pitch line velocity 
on the load capacity of geaurs. The resulting report indicated a pro- 
cedure to determine the so-called dynamic load increment due to 
dynamics of gears in mesh and the error of the gear teeth. Buckinghaun 
presented the dynamic load increment calculation in his text^^^ as 
follows : 

- F+/f“[2F2 - F^] ...(3) 


\diere 


F^ - F[(e/D) + 1] 


F F 
12 

F + F 
1 2 


..( 4 ) 
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= instantaneous load, lb. 

F = average transmitted load (cal ciliated from the horse- 
power to be transmitted cmd considering the force to 
act tangential to the pitch circle) 

F = acceleration load on gear teeth, lb. 

A 

= force required to deform the teeth the amoiant of the 
effective error, lb. 

e = measured error in action (maximum) , in. 

D = displacement of gear tooth under load F, in. 

F^ * force required to accelerate the masses of the gear 
and pinion as rigid bodies, lb. 

The instantaneous load determined by equation (3) should be less 
than the safe allowable load determined from the Lewis Equation. 

Probably, the roost important finding by the ASME Research 
Committee testing program was that most gear failures were not due 
to insufficient bending strength in gear teeth. In many cases teeth 
failed in wear, primarily by progressive pitting. Again, Buckingham 
developed the wear equation which is used today in modified form. 

Tuplin was one of the first to piblish a more refined method of 

[ 2 ] 

det:ennining the dynamic loads in gear teeth. He considered an 
equivalent spring-mass system as shown in Figure 6 that represents 
gears in mesh. He states that passage of a "high" tooth through the 
meshing zone is equivalent to the rapid insertion of a thin wedge 
between loaded teeth of stationary gears and that the model in Figure 
6 represents this condition. The mass M is determined from equivalent 
masses of gears concentrated at the gear pitch cirules. Spring stiff- 
ness K is that of two teeth acting together and is determined from the 
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Figure 6 - Tuplin * s Model 
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static load-deflection relationship. In this model, e and V represent 
the maximum pitch error and rate of insertion of the wedge (average 
pitch line velocity of gears) , respectively. 

Tuplin's res^llts indicate the dynamic load increment to be 


I « 


I ■ 


Ke 

1 + 6.6(t/T)^ 

0.815 X Ke 

[1 + 6.6 (t/T)^]^ 


for t < 0.3 


for t > 0.3 


...(5) 


...( 6 ) 


where 

I > dynamic increment of load (the load above the average 
load), lb. 

K » spring constant of teeth in mesh, lb. /in. 

e » maximum pitch error, in. 

T “ natural period of vibration of equivalent spring-nass 
system, sec. 

t “ time of insertion of wedge, sec. 

It is especially interesting to note that the dynamic increment 
determined by Tuplin's equation has no relation to the average load 
being transxoitted between gears. Also, the equations do not account 
for multiple tooth contact or damping in the system. 

Attia^^^ performed experiments to determine the actual instantaneous 
loaui. He fotand that Buckingham's equation gave high values of dynamic 
increment while Tuplin's equation gave vadues nearer to those measured. 
Also, Attia's measurements illustrated sudden rises and drops in the 
load curves, indicating that gearing errors caused several in?>acts 
throughout engagement rather than smooth load transmission, and that 
the mayiTmim load did not occur at a particxilar phase of engagement. 

From the results he inferred that the simple analyses presented by 
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Buckingham and Tuplin were not adequate to describe the transmitted 
load behavior. 

[4] 

Keswick conducted a more rigorous analyses by incliiding the 
effects of multiple tooth engagement. He also considered the effects 
of heavily and lightly loaded gears. 

Niemann and Rettig^^^ found in their test program that larger 
masses caused higher dynamic loads , but that as the average load 
becazDe larger the effect of larger nuisses bec 2 une unimportant. Tliey 
also found that "very heavily loaded" gears showed no appreciable 
dynamic load increment, %rtiereas in "lig^itly" and "moderately loaded" 
gears dynamic loiui increments of considerable magnitude were observed. 

Harris cwried out a photoelsistic investigation concerning 
dynamic loads of gears. He concluded that when spur gears are isolated 
from external forcing functions, the dynamic load is caused by 

1. EIrror in the velocity ratio measured under the working 
load (gears can only approadi a constant velocity ratio 
under one deflection which depends on the applied load 
and profile modification) . 

2. Pareunetric excitation due to the stiffness variation 
of the teeth. 

3. Nonlinearity caused by tooth separation (backlash) . 

[ 7 ] 

Munro's work in gear dynamics indicated that transients do not 
decay as quickly as previously thought. Hence, he strongly suggests 
that single tooth studies are inadequate, since essential nonrepetitious 
errors are considered. He fovind that after a tooth with error had 
passed through an engagement cycle, subsequent engaging teeth were 
eiffected by the preceding tooth's error* 
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rsi 

Richardson^ cooqpleted an analyses of static load, stress and 
deflection cycles of gear teeth and siobstantiated his results with 
experimental measurements* He then developed a dynamic model to pre^ 
diet the instantaneo\ 2 s load by first considering two gears in mesh ais 
shown in Figure 7. Newton's laws of motion were applied to this 
physical system euid then the system of equations were transformed* The 
model shown in Figure 8 is the result of the transformed equations* 
Assumptions made in order make the problem of determining the instan- 
taneous transmitted loads tractable were: 

!• Input and output toirques remadn constant and the output 
torque is inversely proportional to the velocity ratio. 

2* the total mass M of the model is determined considering 
the equivalent mass of the gears concentrated at the base 
circles of the gears. 

3. Coulomb friction is considered negligible. 

4* The viscous friction force and all other friction are 
considered as a single dancing term. 

5. ^e stiffness of the gear teeth is assumed to be the 
same for all teeth and constant. 

6. Error functions act as forcing fmetions on the system. 

7. Ihe cam moves at the pitch line velocity of the ge£u:s. 

8. Curved ends of the cam result from a "ho load sepeuration 
analysis" as described by Richardson. (The attempt is 

to make the gears engage gradually, rather than abruptly) . 

Two modes of dynamic operation were considered: 

1, Heavily loaded gears where the relative displacement or 
static displacement of gears is greater than the errors 
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Flgvire 8 - Richardson ' s Model 
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involved. 

2. Lightly loaded gears where the errors are much greater 
than the relative displacement, single tooth action 
predominating . 

Richardson's work, for the most part, substantiated each mode of 

operation. His equations are presented in such a manner that the only 

forcing functions acting on the model au:e error functions; however, 

his error function did not aJ.low for the condition of tooth separation. 

It is felt that tooth sepeuration is important in lightly loaded systems, 
roi 

Kastiba analyzed heavily loaded geeirs, and developed a model 
similar to Richardson, but rather than using a cam and an error to 
impress displacement upon the spring, he used a simulated engagement 
error, s(t) . His model, after conversion to a spring-mass system, is 
shown in Figure 9. He also considered a planetary system and used it 
both for analytical and experimental investigation. For the plane teu^ 
system he en^loyed the model shown in Figure 10. The model representing 
the planetary system accounts for the reaction of the system to errors. 
This is the first attempt to consider system effects. Dynamic load 
factors foxmd by Kasuba are generally smaller than those obtained pre- 
viously. He recommended that the actxaal contact ratio under load 
should be used. The entire system of gears should be considered when 
attespting to determine the instcinteuieous load to which teeth are 
subjected. Also, the tooth stiffness should be considered a variable 
accounting for miiltiple tooth contact- 

Bollinger considered tooth stiffness as a trapezoidal 
function. Results of the study correlated very well the experimental 
and analytical work. He found that under different riinning conditions. 


19 


W 



Figure 9-Kasuba*s Single Degree of Freedom Model 
Corresponding to a Compound Differential Gear Train Model 
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the trapezoidal stiffness function for the same pair of gears changed. 

He determined the stiffness function experimentally and used it in his 
analyticail investigations. 

Many investigations have been made to determine the gecu: tooth 
deflections. Baiad and Peterson proposed equations for the deflec- 
tion of gear teeth, considering the tooth as a cantilever beam of 

[121 

variable cross section. Walker deduced a formula for the deflection 
of the gecu: tooth following an experimental, investigation caunried out 
on one tooth fixed to a frame. Both investigations assumed the tooth 
to behave as a cauitilever beam fixed to a rigid wheel. Weber 
considered the actxial shape of the tooth profile in his analysis using 
strain energy techniques. He accounts for normad^ sheair, and bending 
energy in the tooth and a small surrounding area of the gear to which 
the tooth is considered attached. He also considers Hertzian deforma- 
tion assuming tooth profile radii of curvature as equivalent cylinders. 

[141 

Attia expanded Weber's model by including the circumferential 
deformation of the gear rim, auid the deflection of a tooth under the 
effect of loaded neighboring teeth. Further improvements were made by 
Cornell in the treatment of the fillet/foundation deflection. He 
defined the deflection of three different fillet configurations, vdiereais 
previous stxidies assumed a given fillet auigle of 75^. He evaluated the 
resxilting conqpliance anadysis against available test, finite element 
and exact transformation analyses, and foiind that the calculated 
compliance resiilts agreed well with measxirements . Premilhat, Tordion 
and BajOTiet^^^^ determined the elastic compliance through the use of 
appropriate stress functions resulting from the complex transformation 
of the tooth profile. Oheir anadysis resulted in a slightly larger 
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compliance than Weber's analysis. Chabert, Tran and Matis^^^^ made 

yet a different evalxiation of stresses and deflection of sp\ir gear 

teeth by finite element methods. Formulas were developed for a 

simple calculation of the maximum stresses and compared the results 

with ISO and AGMA standard formulas for the strength of gear teeth. 

[181 

Cornell and Westervelt ^ developed a time history, closed form 

solution of a dynamic model of spur gear systems. Their emalysis 

determines the dynamic response of the geeu: system and the associated 

tooth loads and stresses. The dynamic model is based on Richardson's 

cam model but treats the teeth as a variable spring. Included in 

the analysis are the effects of the nonlinearity of the tooth pair 

stiffness dxiring mesh, the tooth errors, and the profile modifications. 

The anailysis showed that tooth profile modification, system inertia 

and damping, and system critical speeds can affect the dynaunic gear 

tooth loads and stresses significcuitly. 

[191 

Kasuba and Evans concluded that the gear mesh stiffness in 
engagement is probably the key element in the analysis of gear train 
dynamics. Also, the gear mesh stiffness and contact ratio are affected 
by many factors such as the transmitted loads, load sharing, gear tooth 
errors, profile modifications, gear tooth deflections, and position of 
contacting points. 

By introducing these aspects, the calculated gear mesh stiffness 

can be defined as being a variable - variable mesh stiffness (WMS) . 

The WMS model is an improvement over the previous periodic varying 

[191 

mesh stiffness which the authors called fixed - variable mesh 
stiffness (FVMS) . 

They developed a large scale digitized gear model including the 
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WMS method which investigates in one uninterrupted sequence the static 
and dynamic conditions • 

An iterative procedure wais \ised to calculate the WMS by solving 
the statically indeterminate problem of multi-pair contacts , cdianges in 
contact ratio, and mesh deflections. The developed method can be used 
to ancd.yze both normal and high contact ratio gearing. 

The associated computer program package calculates the WMS, the 
static and dynamic loads, and variations in transmission ratios, sliding 
velocities, and the maocimum contact pressures acting on the gear teeth 
as they move through the contact zone. 

Their findings for typical single stage external spur gear 
systems are: 

1. The geaurs and the adjacent drive and load systems can be 
matched for optimum performance in terms of minimum 
allowable dynamic loads for a wide range of operating 
speeds. 

2. Torsionally flexible design bf geair bodies/hiabs/rims 
offers an excellent means of absorbing or minimizing the 
geometrical errors in mesh. 

3. The geau: mesh stiffness and its distribution axe 
significantly affected by the transmitted loads and tooth 
profile imperfections. 

4. The dynamic factors can be decreased by increasing the 
damping and/or contact ratio. Local damping appears to 
be the most efficient means for decreasing the dynamic 
load factors. 

5. The high contact ratio (HCR) geairing has lower dynamic 
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loads and peak Hertz stresses than the normal contact 
ratio (NCR) gearing. 

The studies by Kasuba and Evains provide one of the most detailed and 
advanced models availedsle at this time. 
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2.2 INTERNAL SPUR GEARS 

Buckingham indicated that there are almost an infinite number 
of forms v^ich can be used as gear tooth profiles. However, the most 
common profile for transmission of power is the involute form. He 
developed the kinematic formulations under which involute gears trans- 
mit uniform rotary inotion. The constant velocity action between such 
teeth is called conjugate gear tooth action. Buckingham's formulations 
of the conjugate gear-tooth action and interference prediction for the 
internal spur gecu: are readily applicable for present day use. 

Dudley further pointed out that a comparison between the 
internal and external spur gear sets assuming the same number of teeth 
produced the following advemtages of the internal spur gear set: 

1. Greater length of action. 

2. Relative sliding of the teeth at the start and end 
aire less. 

3. Center distance is smaller and thus leads to a more 
compact arrangement. 

4. Contact area is larger because of the mating of a 
concave and convex surface. This increased contact 
curea results in larger resistance to pitting and 
wear. Also, the distribution of the load among more 
teeth decreases the intensity of the stress. 

An additional advaintage is derived because the tooth strength of am 
internal gear is greater than that of an equivalent external gear. 

Dudley acknowledges that there is no AGMA standard covering 
internal gears. He recommends that the methods of design for the 
external spur gear may be applied to the internal spur gear . In light 
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of the obvious advantages of the internal set, it is understandable 
that the external spur gear design methods would lead to a satisfactory 
internal gear design. Incidentally , the available methods for external 
gears that Dudley refers to ajre primarily based on static analysis of 
the system. In this study the internal spur gear set will be analyzed 
from the dynamical point of view. 

Dudley also lists several disadvantages to the internal gear set. 
These disadvantages can be removed by introducing more accurate and, 
consequently, more costly manufacturing operations. 

The first of these disadvantages is tip interference. In this 
type of interference, the external gear cannot be assembled radiailiy 
with the internal gear. Only axial assearibly is possible. If a shaper 
cutter having a n\mber of teeth equal to or greater than the external 
gear is used to cut the internal gear, then it will cut its way into 
mesh but in so doing will remove some material from the flanks of a 
few of the teeth that should have been left in place for good tooth 
operation. This cutting action is also known as ’'trimming". Such 
teeth will have poor contact and will tend to be noisy. 

The second problem is sometimes known as "fouling". In this 
case the internal gear teeth interfere with the flanks of the external 
toothed gear if there is too small a difference in numbers of teeth 
between the external and internal gears - 

[ 21 ] 

The third problem is in the manufacture of internal gears. 

The necessity for the generating tool to work within the gear body 
restricts cutter and machine dimensions, w’hich in turn limits tool 
accuracy, rigidity and the resultant precision. Finish grinding of the 
teerh is especially difficult because of the large size of the grinding 
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wheel. Generally, internal gear sets are not used for precision geau: 
trains unless the applications must utilize their unique internal fea- 
tures as an advantage. Examples of such applications are indicated in 
Section 1.1, "General Remarks". 

Clearly, any detailed investigation of the ISG drives, such as 
this thesis work, must consider the previously mentioned findings and 
recommendations by Buckingham and Dudley. However, the above findings 
must be used in conjunction with the modern thinking of gear behavior 
as was indicated in the literature search for external gears. This 
modern thinking must include consideration of the elastic deformation 
of the internal gear teeth and their supporting ring structure, and 
dynamic analysis of the system. 

To date, the number of investigations related to the elastic analy 

sis of the internal spur gear teeth and the supporting ring are limited 

[ 22 ] 

Karas was the first to evaluate the deflection of the internal 
spur gear tooth due to bending, shearing and Hertzian contact deforma- 
tion. He assumed a trapezoidal shape for the tooth profile and a rigid 

support at the root. He did not consider any deformation of the sup- 

[23] 

porting ring, Ishikawa regarded the root of the gear tooth as a 

semi-infinite body, and then calculated the deflection due to the 

tooth rotation at the root of the gear tooth. He also did not consider 

[24 ] 

the deformation of the ring. Hidaka superimposed the above four 
deflections into one final tooth deflection. He then compared this 
deflection against a finite element representation of the tooth auid 
portions of the ring. The boundaries of the ring portions were then 
fixed against translation or rotation. He concluded that the results 


of his four deflection relations asid the finite element analyses were 



relations for deflection due 


, , [25 

L uy o2.j;iKe vx c:n . 


similar. Thus, he decided to use the four 

to bending, shear, Hertzian contact stress and rotation at the root in 
the analysis of planetary spur gear systems. 

The fir 

He replaced the ring gear with a perfect ring having an egui'/alent stiff- 
ness. The equi'^/alent thickness was expressed as a function of module/ 
diametral pitch, whole depth of the tooth, backlash and the number of 
teeth . 

[24] 

Hidaka later compared Sinkevich’s relations agacrmt a fcrxite 

element representation of the ring. He treated the deflaction as a 
plane stress problemi usxng different finite element mesh si::::es and 
different thickness. Having arrived at an optimum size and thick- 

ness he then modified Sinkevich’s deflection relations based on the 
equivalent thickness concept. However, Hidaka 's final L(d.aticns for 
ring deflection are not applicable to the present ISG ri.nc gear deflec- 
tion investigations due to the following reasons: 

1. His investigation is for a planetary gear system xn 
which the loading is symmetric around the circum- 
ference. The ISG drive has a single point loading. 

2. He assumed thin ring relations (i.e., thickness over 
radi.us ratio 1/10) . In many cases the ring gear of 
Che ISG drive is a thick ring. 

Hidaka also points out that the finite element method of solution for 
the deflection of gear teeth requires a finely meshed model . This 
approach can quickly exhaust the capacity of the computer. 

Because of the limited available information for deflections of 
the ISG drives, the author of this thesis decided to utilize the 
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appliccdile and proven methods for determining the stiffness of external 
gearing systems to the ISG drive. The work by Weber and Attia^^^^ 
on external gear teeth and hubs will be adapted to the internal tooth 
profile and then will be used as a basis for comparison of resvilts. 
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CHAPTER III 

ANALTOICAL INVESTIGATION 
3,1 PROBLEM FORMULATION 

A comprehensive analysis of the dynamic loading in interned spur 
gear (ISG) drives presents a difficult task even under ideal geometry 
conditions becaiise of the continuously changing interactions as the 
gears move through the meshing cycle* Further complications arise 
when the gear tooth deflections* backlash* profile errors and multiple 
gear tooth pair loading are introduced. The complexity of the meshing 
process can be illustrated by following the action of the teeth through 
one complete mesh cycle. Figure 11 shows the two gears at the start 
of the cycle where the driver engages the driven tooth at its tip. 

Onder load the teeth deform to a noninvolute shape which changes the 
line of action and* thus, the loading between the teeth. Other in-plane 
deformations, both teuigential and radial , taUce place in the tooth support 
structure, the adjacent teeth and the remaining drive train. Figure 12 
depicts this deformed condition. If the deflection continues, then 
contact will be made between another tooth pair and an indeterminate 
load sharing condition is entered. Backlash or a tooth that is too 
thin increases the deflection slightly but In general decreases the 
clumce of multi-tooth contact. A tooth that is too thick leads to 
premature engagement cuid jamming between teeth. A pit in the profile 
can cause sudden disengagement and subseqxjent clashing between teeth. 
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As the driver pushes through the mesh cycle, the loading changes from 
root to tip contact and finally disengagement. The unloaded teeth will 
regain almost all of their original shape immediately, and then com- 
pletely as additional teeth become unloaded. A typical combined 
stiffness pattern for error free teeth is shown in Figure 13. This 
periodically repeating pattern will be distorted by identical profile 
errors in the teeth. The nonlinear dynamic process leads to instanta- 
neous load fluctuations in the teeth even in the presence of constant 
external load conditions. Also, the magnitude of the load and the 
fluctiiations are influenced by the damping effect of the lubricamt, 
and the proximity of component natural frequencies with any of the 
forcing frequencies. Figure 22 shows a practical internal gear drive 
system model used in this study. 

Having thus identified the physical problems it can be stated here 
that adequate mathematical tools are available for their solution. 
However, the capacity of the present large scale digital coir^uters and 
the scope of this investigation are such that certain limitations must 
be imposed on the treatment of the problem. These limitations or 
assimptions are treated next. 




Figure 11 - Internal -Extemail Gears at Start of 
Mesh Cycle Under No-Load 
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Figxare 12 - Internal-External Gears Under Load 
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Figvire 13 - Combined Stlffiiess Pattern of External-Internal Spur Gears 
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3.2 ASSUMPTIONS 

The following sinqplifying assumptions and conventions are used to 
make the 2 Uialysis manageable and still realistic: 

1. The dynamic process is defined in the rotating plane 
of the gears. The torsionad stiffness of the shafts 
cuid gears in engagement, and their masses are acting 
in the plane. This assumption is considered realistic 
because of the symmetry of the rotating axes. Also, the 
out-of-plane twisting and misalignment are prevented by 
proper design and careful assembly procedures- 

2. Vaiaping due to Ixabrication of gears and bearings are 
expressed as constant damping coefficients. Their 
effect on the load dynamics is investigated by para- 
metric studies. 

3. The dynamic process is investigated through several 
complete engagement or mesh cycles. 

4. The differential equations of motion are expressed 
along the instantaneous rather than the theoretical 
line of action of the teeth which permits evaluation 
of noninvolute action. 

The deformations of the tooth support structure and 
shafting are determined from equations which were 
developed in solid mechanics for simple shapes. 

For more complicated shapes the deformations have 
to be determined experimentally or by finite element 
techniq\aes. The con^uter program is structured so 
that the deflection values from experiment or finite 
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element analysis can be entered into the analysis 
process by means of data sets. 

6 . Multi -tooth contact is determined by analyzing five 
gear tooth pairs. The central or middle tooth is used 
to establish the instantaneous position of the teeth 
and monitor their progress through a complete mesh 
cycle . 

7. The presence of backlash may lead to tooth separation 
under dynamic conditions which must be accounted for 
in the emalysis methodology. 


37 


3.3 METHOD OF SOLUTION 

The purpose of this investigation is to develop am analytical 
method for determining the static amd dynamic behavior of the ISG 
drive. In order to obtain dynamic information, it is first necessary 
to get the static supporting data. For this reason, it has been con- 
venient to divide the investigation into static and dynamic sections. 

In each section, the appropriate emalytical model is comprised of 
relations available from references in gearing, strength of material, 
mathematics, vibration analysis, and the publications on deflections 
by Weber, Attia, etc. In the interest of timely solutions, an 
attempt was made to solve for the required information directly. Where 
this was not possible, iterative search techniques and numerical solu- 
tions, along with suitable acceptsmce criteria, were substituted. 

The developed analytical methods were combined in a sequence of 
digital computer programs which can be used on a large scale computer 
like the IBM 370/158 at Cleveland State University (CSU) . For parametric 
studies, the program can be used more efficiently in three parts (modules) . 
The first module determines the static information and stores it on a 
tape. The second module uses the static data to initiate the dynamic 
solution and then solves for the dynamic information. The third module 
calciilates the static or dynamic tooth bending stress at the critical 
fillet location. In a similar fashion, the program can be fiirther sub- 
divided for incorporation on a mini-computer like the Hewlett-Packard 1000. 

The next three sections of this report present a detailed develop- 
ment of the static and dynamic analysis procedures as well as a s\ 2 mmary 
description of the computer program. For reference purposes. Appendix A 
cmd B contain all of the standard eqiiations which were utilized in this 
investigation. Appendix C contains the computer program listings. 
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3.4 STATIC ANALYSIS 

The task of the Static Analysis Procedure is twofold. First , it 
must provide all of the supporting and final information needed from 
a static analysis of a gear system. The emalysis procedxire must be 
structured so that the desired information is obtained efficiently. 
Some of the structxaral requirements and the needed information can be 
identified as follows: 

1. A suitable nomenclatxire for documentation and 
computer \ise. 

2. SuiteQ)le local and global coordinate systems. 

3. External and internal gear tooth profiles. 

4. Contact points between gear tooth pairs. 

5. Line of action. 

6. Contact ratio. 

7. Interference conditions. 

8. Deflection and stiffness of the teeth and their 
supporting structure. 

9. Load shairing among neighboring teeth. 

10. Sliding action between mating teeth. 

11. Static load per tooth pair. 

Second, the static analysis must file this information for use 
in the dynamic and stress analyses, and for printing of selected 
portions of this information. 

3.4.1 Nomenclature 

The nomenclature for the static and dynamic analysis has been 
selected from symbols used in gearing, strength of materials, mathe- 
matics and publications by Weber, Attia, etc. When the required 
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symbols were not available from these sources, special symbols were 
introduced to describe the particular parameters in short form. 

3.4.2 Local and Global Coordinate Systems 

Three Cartesian coordinate systems are employed in the static 
analysis. The first is a local system using the symbols X and Y. It 
has its origin at the root of each tooth. The Y-axis coincides with 
the tooth centerline. In all, there are ten such X-Y local coordinate 
systems to account for five gecu: tooth pairs under investigation in 
each gear. Treuis formation from one tooth coordinate in a given gear to 
another tooth is readily possible because of the fixed geometric rela- 
tion between the teeth. These local tooth coordinate systems are used 
to define the discrete tooth profile locations and appropriate deflec- 
tions of the teeth. 

The coordinates of the second system are laO^eled W and Z. These 
systems are local to each gear and rotate with the gear . Each W— Z 
coordinate system is parallel to its respective X-Y system. There are 
ten W-Z coordinate systems also. The Z-axis coincides with the tooth 
centerlines. 

The third system is global and fixed at the center of the internal 
gear. This system is identified as the U-V system. The arrangement 
of the three coordinate systems is shown in Figure 14 . The transfor- 
mations between the coordinate systems for each gear pair are: 

W1(I) * XI (I) 

W2(I) * -X2(l) 

Z1(I) « RROl + Y1(I) 


Z2(I) * RR02 - Y2(I) 
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U1(I) = W1(I) X SIN[PSI1TP(I)] + Z1(I) X C0S[PSI1TP(I)] 
V1(I) = -W1(I) X C0S[PSI1TP(I)] + Z1(I) X SIN[PSI1TP(D] 
+ C 

U2(I) = -“{Z2(D X SIN[PSI2TP(I) - 0*5ir] - W2(I) x 
C0S[PSI2TP(I) - O.Stt]} 

V2(I) = Z2(I) X COS[PSI2TP(I) - O.Stt] + W2(I) x 

SIN[PSI2TP(I) - O.Stt] ...(7) 

where 

I » gear pair 1 through 5 

1 = external gear 

2 * internal geatr 


3.4.3 External and Internal Gear Tooth Profiles 

Development of the involute profile of a tooth follows well 
known geometric relations which are shown in Appendix A for convenient 
reference. Also shown in Appendix A are construction of the involute 
profile, and the respective external and internal local and global 
coordinates of the computer program. 

In actual practice, deviations from the theoretical involute 
profile are introduced because of manufacturing tolerances, errors 
and by intent. A fourth potential source of modification is due to 
damage during operation. The intentional modifications are introduced 
to overcome the detrimental effects of profile deviations and toler- 
einces in the remaining system. As a rule, gears are fcO^ricated with 
backlash on their teeth. Additional modifications are made to the 
tip to improve engagement between mating teeth, and to the root to 
avoid interference. Damage during operation exhibits itself as local 
burnishing, pits or spalling. The condition of a tooth profile can be 
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determined by means of a gear checker. This machine follows the tooth 
s\irface as the gear is rolled of its base circle. A moving pen and 
cheirt instrument then draws the profile as a function of the roll angle. 
A true involute produces a straight line and cuiy deviation from this 
line is indicative of the degree of modifications and/or faults. 

Thus, it is customary in gear analysis to define the profile 
modifications ctnd errors by means of a profile chart. Figure 15 illus- 
trates the relationship of the profile and involute chart for the 
external spur gear tooth. As indicated before, the profile modifica- 
tion of the tip deviates from the straight line of the true involute. 

A simileur chart holds for the internal spur gear tooth. Figrure 16 
depicts samples of possible tooth profile modifications which can be 
used alone or in combination to simulate a large number of practiceO. 
cases. Various profile modifications for the external and internal 
tooth are shown in Figxires 17 and 18 respectively. 

Analyticailly, the profile modifications can be expressed by meauis 
of a prodiict consisting of the maximum cunplitude of the modification 
and an appropriate shape function. This shape function uses the roll 
auigle as the independent variable which can be structured so that the 
straight-line, parabolic, sinusoidal and pit like modifications of 
Figxire 16 are faithfxilly represented. In this manner, the surface 
faults of a particular profile are combined with the true involute 
profile, as developed in Appendix A, to simulate the actual tooth 
shape. The X-Y coordinates of such a tooth profile from the tip to 
the beginning of the fillet are defined by the following general 
expression: 


I 



TIP AND ROOT HODIFICATIONS 
(STRAIGHT LINE OR PARABOLIC) 
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Figure 16- Sample Simulated Geau: Tooth Profile Charts 
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External Tooth (sijbscript 1 omitted) 


X = R SINB ± 


STTM 

cose 


RA - RAM ^ PATH 
^ RATM ^ ’ cose 


[1 - { 


RAT - 


RATM 




* 55 se 5 m » ’’’^1 - 

^ STEM RA - RAN ^ PABM RA - RABI ^ 

cose ^ RABM ^ “ cose " RABM ^ ^ 


Y = R C0S6 “ RRO 


...( 8 ) 


Internal Tooth (subscript 2 omitted) 


X 


-{R SINS 


^ SiTM 

“ cose 


RA - RAM . 
RATM ■ 


^ PATM , _ RAT - RA ^ 

cose ” ^ RATM 


+ PER 

“ cose 


SIN [7T(RAM - RA) 


CYC 

RATIP 


+ PAP] 


DEEP 


^ STEM , RAN - RA ^ PAEM , RAEI - RA ^ 

* cose ^ raem ^ “ cose ^ ~ ^ rabm ^ 


where 


STTM » magnitude of straight line modification at tip 
RA > roll angle 

RAM = roll angle at end of modification at tip 
RATO = length of tip modification in degrees of roll 
PATM = magnitude of parabolic modification at tip 
PER * maximum profile error 
CYC = number of sinusoidal error cycles 
RATIP= roll angle from end of modification to roll 
angle at the pitch point 

PAP = angle from start of sinusoidal error to start 
of involute 
DEEP * depth of pit 

STEM = magnitude of straight line modification at bottom 
PAEM = magnitude of parabolic modification at bottom 
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RABM = length of root modification in degrees of roll 
RABI * roll amgle at the bottom of involute 
RAN ^ roll angle at end of modification at bottom 
and the remaining paraimeters are defined in Appendix A, In the computer 
program, the tooth profile is represented by a finite number of points- 
The spacing between points is selected so that the segment between two 
points can be represented by a streuLght line. Depending on the size of 
the tooth, either one or two hundred points are used. 

Thus by specifying the appropriate parameters euiy profile config- 
laration of Figures 17 and 18 cam be represented by equations (8) amd (9) 
in digital form. The development of the fillet coordinates is shown in 
Appendix A. 

3.4.4 Contact Points Between Gear Tooth Pairs 

The location of the contacting gear teeth and the number of 
contacting gear tooth pairs caumot be determined directly by analytical 
means because of the possible presence of tooth modifications and the 
deformation of the loaded teeth. However, it is possible to determine 
the line of action, and initial and final points of contact of teeth 
with true involute profiles. Figure 19 shows the arramgement for an 
external-internal gear pair. The theoretical initial and final points 
of contact as shown in Figure 19 are used in this investigation as the 
starting points for a two step iterative seairch of the actual contact 
points as the gears rotate through a complete meshing cycle. The 
search for this actual contact considers first the unloaded (rigid 
body) motion of the gears and then repeats this search again while 


the gears rotate under load: 



LINE OF 
ACTION 
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Step One - Rigid Body Motion 

1. Transfer the local X-Y coordinates of the digitized tooth 
profiles to the rotating W-Z and global 0-V coordinates. 
This transfer fixes the tooth (assumed to be the third of 
five) unto the respective gear and establishes the geome- 
tric arrangement between the two gears. 

2. Determine the theoretical initial amd final points of 
contact from the standard gear relations of Appendix A 
wd the geometric arrangement of the two gears as shown 
in Figure 19. 

3. The theoretical initial and final position of the center 
of the tooth are determined. The angular arc between the 
two angles is divided into forty-nine equal intervals. 
Since the true and actual positions of contact for modi- 
fied involute profiles may differ, it is necessary to 
search for these positions. This search is accomplished 
by rotating each gear several interv 2 d.s counter clockwise 
and thoi the proximity of each profile point of gear 1 
and 2 is conpared. When the proximity of any profile 
points of gear 1 and 2 fall within the specified accep- 
tance criterion of 1.0025| mm then actual contact is con- 
sidered to be made. This search is performed at the 
initial and final point of contact. The angular arc is 
again determined and divided into forty-nine intervals. 
This new interval is used as the angle with \diich the 
gears will be rotated in forty-nine increments resulting 
in one conplete mesh cycle. When the proximity of any 
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profile point is not acceptable, th&i the two gears are 
rotated clockwise a fraction of an interval and a new 
comparison is made until contact is found. For unloaded 
teeth with true involute profiles, the contact occurs 
near the theoretical point of contact. For modified 
profiles, this search may continue beyond the original 
theoretical point of contact. In Figure 20, the distance 
(L) represents the closest distauice between the 
tip of the interned, gear tooth and end of the involute 
profile of the external gear tooth. The point is 
determined from the sxirro unding profile points and by 
using the relation between similar right triangles: 


U 


11 


V2(L) - V^(J) 

v^Tj + 1) - v^(j) 


[U^(J + D - U^(J)] + U^(J) 


... ( 10 ) 


The final contact position is determined in a similar 
fashion except here the tip of the external gear tooth 
and the end of the involute profile of the internal gear 
tooth are in contact. 

4. The initial and fined, tooth center position of tooth 
number one is determined by subtracting two circ;ilar 
pitches from the angular position of the third gear 
tooth center. Incrementing in integrad. values of a 
circular pitch results in initial and final positions 
of teeth two through five. At this point the unloaded 
initial and final positions of all five teeth of the 
internal and externad. gear are established, and tooth 
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pair number three is contacting at the beginning of its 
mesh cycle. 

5. A complete contact search for all five tooth pairs is 
made at fifty positions from initial through final con- 
tact of tooth no. 3. Actual contact at every mesh point 
for all five teeth and the contact position for tooth 
number three and any other teeth is recorded. At this 
point single or multiple tooth contact for all fifty mesh 
positions is determined. In case of multiple tooth con- 
tact the load will be shared between the contacting teeth. 
Each tooth carries a fraction of the total load propor- 
tional to its own tooth sti-ffiiess and the combined stiff- 
ness of all the contacting teeth. For the i^^ contacting 
tooth the shared load Q(k)^ is 
KP(k) . 

Q(k)i - ^ (Q^.) ...(11) 

i 

where KP(k)^, KG^ and are the stiffness of the i^^ 
tooth, total mesh stiffness and total load at that mesh 
position respectively. Details of the methods for tooth 
deflection emd stiffness calculations are discussed in 
Section 3.4.6. From these methods and knowledge of the 
shared load on each tooth it is now possible to deter- 
mine the deflection of each tooth profile point. This 
completes the first search for contact assuming rigid 
body motion of the two gears. 

Search for Contact of the Loaded Gears 


1. The gear tooth deflections can be considered as another 
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form of tooth modification causing premature engagement 
and delayed disengagement. Points A* and B* of Figure 
21 demonstrate this action. Thus, by adding the tooth 
deflection, the whole procedure of the first search 
routine must be repeated. At the end of the procedure 
the actual contact positions, mesh stiffness and static 
loauiing for fifty positions of the total meshing arc 
are determined. Also calculated are the deflections 
of each profile point for all five tooth pairs at each 
of fifty mesh angles. 

3.4.5 Line of Action, Contact Ratio amd Interference Conditions 

The line of action of an unloaded involute tooth pair is tangent 
to the base circle of gears 1 and 2, and its inclination is represented 
by the theoreticed pressure angle, (see Figtire 19). Under load the 
mating teeth deflect and the instantaneous line of action is no longer 
tangent to both circles. The line of action will change with changes 
in load, speed aind position. Thus, we are concerned with the theoretical 
and instantaneous lines of action (see Figure 21). Other parameters 
that Ccm be thought of as chainging at any instant are pressure angles, 
base radii, pitch radii, transmission ratios, etc. The net effect of 
the change in line of action is to reduce the transmission ratio, TR, 
or the ability to transmit the torque. 

As mentioned before, the circumferential deflection of the teeth 
causes premature engagement and delayed disengagement. This condition 
is beneficiea to the action of the gears since it increases the arc of 
contact between the mating teeth. In Figure 21 the angular arcs A-B 
and A*-B* represent the theoretical and instantaneous contact arcs. 
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Thus the loaded contact ratio, CR, for an error-less gear pair can be 
approximated as 



^ere 

p* m circular pitch under load 

The deflection of the gear tooth under load enters into another 
consideration, namely, the possibility of interference. Formulae for 
avoiding the various interference conditions of external-internal in- 
volute geeurs have been derived from unloaded gears in terms of limit 
radii or angles* . For the loaded condition the deflection of the teeth 
may increase the potential for interference. Thus, it is necessary to 
conduct two interference calculations using the ideal and instantaneous 
gear geometries* 

Of course, any deviation from the involute action will increase 
the amoisit of sliding action between the teeth* 

3.4.6 Deflections and Stiffness of the Teeth and Their Supporting 
Structure 

When the teeth of a pinion and gear come into contact, trans- 
mission of load caxises deflections in the gearing system as indicated 
in Figure 12* For the external -internal gear combination the overall 
deformation may be sufficiently large to influence the mesh character- 
istics and, thus, affect the static and dynamic behavior of the gear 
system. The actual deflection may be considered to consist of the 
susnation of several deflection oomponents starting with the Hertzian 
deformation at the contact point and extending into the foundation of 
the system* For the dynamic analysis it is necessary to consider 
inertial and damping effects as well which can best be treated by 
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siabdividing the gearing system into smaller parts of masses, dampers and 
springs. In the context of this system siibdivision, the deflections 
pertaining to the contacting gear teeth are radial and circtimf erential • 
The radial component is due to the radial deflection of the support bear- 
ings and shafts and the internal gear ring. The radial deflection of 
the pinion is negligible due to its rigid construction. Expressions for 

the radial deflection of bearings and shafts are readily available from 
T261 T271 

standard textbooks. A listing of deflection ecjuations of bear- 

ings suitaQ^le for the ISG drives is given in Appendix B. The radial 
deflection of the ring gear is not as readily available because of its 
single point loading and complex geometry (see Figure 3). In general, 
experimental means or three-dimensional finite element analysis must be 
resorted to for an accurate assessment of the deflection. However, 
these measures may not be needed for the ISG drive because it is usucilly 
designed to minimize radial deflection. In Figure 4, the support beauring 
for the ring has been placed directly over the load preventing radial 
deformation of the ring. In Figure 3, the ring thickness over the inter- 
nal geeu: teeth is heavy causing nearly rigid body loading of the outer 
cylindrical section. For this configxiration a worst case condition can 
be assxixned to occur as outlined in Appendix B-2 . The worst case solu- 
tion of Appendix B-2 is then used to evaluate the performance of the 
gears. Because of the rigid design, the radial deflections are held 
to the saioe order of magnitude as the circumferential deflection. 
Nevertheless, the radial deflection causes radiaJ. movement of the gears 
with a resulting reduction in contact length of the gears. 

For a given design, selection of the appropriate bearing deflec- 
tion equation from Appendix B-1 and bracketing of the ring deflection 
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leads to a combined radial deflection as follows 


6—6 +6 +6 t 6 

R REB RIB RES RR 


...( 13 ) 


vdiere 

6 = deflection of external gear bearing 

REB 

6 = deflection of internal gear beairing 

RIB 

6 = deflection of external gear shaft 

RES 

6 * deflection of ring 

RR 

The deflection along the gear circxjmference is due to Hertzian 
deformation, bending, shear, con^ression and rotation of the teeth at 
their root, and due to torsion of the pinion and qeax ring. 

The coiriDined external -internal gear tooth pair deflections can 
be expressed in the following form: 


6{k)^ - ^ 


...(14) 


\diere 


6. (k) . ■ deflection of the k tooth of the extezneil 
E 1 

geair at mesh arc position i 

6i(k)i ■ deflection of the k^^ tooth of the internal 

gear at mesh arc position i 

6 (k) . » localized Hertzian deformation at the point 
H i 

of contact 

For the contacting pairs, the gear tooth deflections suid 

6j(k)^ incorporate a number of constituent deflections. For the 
extematl gear: 

■ ‘me*"’! * *Se"i * V"’! * ’ ‘rE^’i 

. . ■ ( 15 ) 

In equation (15) , 


6 * gear tooth deflection due to bending moment 

ME 
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» gear tooth deflection due to shear force 

* 9 ^^ tooth deflection due to normal force 

^BE * tooth deflection due to deformation of 

the surrounding hub area (rocking action) 

^RE * 5®®^ tooth deflection due to torsion of the 

rim or hub (circuxnferentiad. deformation of hub) 

For the internal gear, the deflection for the tooth pair at 
mesh arc position i is 

■ 'hi*’! * * W">i * «BlWi * 

where 

« geau: tooth deflection due to bending moment 

* geau: tooth deflection due to shear force 

* gear tooth deflection dxae to normal force 

* geair tooth deflection due to deformation of 
the surro landing ring area (rocking action) 

* gear tooth deflection due to torsion of the 
supporting ring (circumferential deformation 
of the ring) 

Expressions for the constituent deflections of the external geair 
have been derived by strength of material techniqioes 

In this investigation, the same methods have been applied to the 
internal gear. A detailed account of aJ.1 the circumferential deflections 
in the ISG drive is shown in Appendix B--3 

The circumferential deformations of the hub and ring affect the 
deflections on clLI teeth whether they are loaded or not. Thus, if the 
rigid body contact search of Section 3.4.4 does not find contact between 
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two teeth, it is possible that the two unloaded teeth would be declared 
in contact in the second contact search as a result of the attendant 
circumferential deformation. In this case, the final load sharing and 
deflections will be recalculated on the basis of this additional con- 


tacting tcx>th pair. 

3oth ‘the radial and circumferential deflections affect the mesh 
stiffness and contact ratio diaracteristics. The radial deflection 
primarily affects the contact ratio and to a smaller degree the stiff- 
ness value. The circumferential deflection has a much more significant 

on the mesh stiffiiess. These effects are shown in Figures 26 & 29 

Thus, for any mash arc position i, the calculated gear tooth 
pair stiffiaess KP(k)j^, mesh stiffness KG^^, and load faring incorporate 
the effects due to manufactured profile errors, profile modifications, 
and radial and circumferential deflections by means of the iterated 
nmerlcal solutions of eguations (13) through (18) . 

The individual gear tooth pair stiffiaess can be expressed as 


KP(k) 


QOn)^ 

i ” «(^ 


...(17) 


If the effective errors prevent contact, KP(lc)^ * 0. 

The sum of geeu: tooth pair stiffiiesses for all pairs in contact at 
position i represents the variable mesh stiffness KGP 


K 

KGP. - Z KPOc) . 
^1 ^ 


.(18) 


The load carried by each of the pairs moving through the mesh arc in 
the static mode can be determined as 
KP(lc), 


Q(k) 


KGP. 


(P) 


...(19) 
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where P is the total normal static load carried by the gears at any 
mesh position in the static mode 
K 

P = 2 : Q(k) ...(20) 

1 

The static analysis thus described determines the variable mesh 
stiffness (KGP) , transmission ratios (TR) / and the contact position 
vectors (RCPl/ RCP2 , RCCPl, as shown in Figure 21) for subsequent 
dyneunic ceilculations. 
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3.5 DYNAMIC ANALYSIS 

Itie gear train sho\m in Figure 22 was modelled for the dynamic 
solution. This gear tradji is found in practical applications like t\ir- 
bine driven pumps, motor driven tank turrets or wind tiurbines as in 
Figure 4. The model consists of input and output devices, the external- 
internal gear transmission and interconnecting shafts and bearings. 

The emalysis considers constant input and fluctuating output torque, 
damping in shafts, geaurs and bearings, backlash, noninvolute action 
caused by deflections and tooth modifications, and loss of contact be- 
tween gear teeth. The coordinate system used in the static analysis 
is used also in the dynamic analysis. The instantaneous parameters 
^ich were determined for the fifty mesh arc positions in the static 
analysis will be combined with the equations of motion for the dynamical 
solution of the system. 

For the model of Figure 22, the differential equations of notion 


cam be given in the following form: 

«D * ‘=B0 «D t ' ®2> 


* - “l> ■ ■'c 
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"'l «L * =BS \ * =BL * =LS<®L ' * \s<«L ‘ ®2> 


= - T_ X TR' = - T_ (TR') 
D Xi 


... ( 24 ) 


Importauitly, the equations of motion are based on the instantemeous 
rather than the theoretical line of action. In equation (24) , the load 
torqxje is written as a function of the instantaneous transmission ratio 
TR* . The bracketed terms in equations (22) and (2 3) represent the dy- 
namic gear mesh force ^^ich is dependent on the dynamic displacements 
of the engaged gears, geau: mesh stiffness and damping in the mesh. 

The mesh stiffness, KGP^ in equations (22) and (23) , represents 
the combined effects of gear tooth profile errors and modifications, 
radial and circumferential deflections of the gear teeth, sharing of 
the load between teeth, height of engagement, and the angular position 
in the gear mesh cycle. Representative mesh stiffness cycles are 
shown in Figures 26 and 27. For a constant input torque, the variable 
mesh stiffness and the changes in transmission ratio due to noninvolute 
action represent major sources of dynamic excitation. The output torque, 
T_ , is a function of the input torque, the instantaneous transmission 

ii 

ratio and losses in the system. 

Diiring operation of the system in Figure 22, the dynamic excitation 
sources C2in create sitioations during which momentary disengagement of 
the mating gear teeth can occur. The information vdiether separation 
takes place can be obtained by reviewing the equations of motion. 

The term (RBCl 0^^ - RBC2* 0^) represents the relative dynamic 
displacement of gear 1 and 2. Considering that gear 1 is the driving 
gear, the following situations can occur: 



INTERNAL GEAR 


64 



Figure 22 - Internal -External Gear Train Used in the Dynamic Analysis 



65 


When RBCl 0^ > RBC2* 0^ 

we have normal operation of the gear system and the dynamic mesh 
force is defined by 

(gOT)^ * CGP^(RBC1 0j^ - RBC2* 0^) 

+ KGP^(RBC1 0^ - RBC2' 0^) -..(25) 

If RBC2 0^ ^ RBCl* 0^ and RBC2* 0^ - RBCl 0^ < BGM 

then the gears will sepaurate and oontact between the gears will be 
lost. For this case, 

(QDT)^ - 0 ... (26) 

If RBC2* ©2 > RBCl 0^ and (RBC2* 0^ - RBCl 0^) > BGM 

(QDT)^ * CG^(RBC1 6^ - RBC2* 

+ KGP^KRBCl 0^ - RBC2* - BGM] ...(27) 

in this case, gear 2 will hit gear 1 on the backside. 

Also when KGP^ * 0; (QDT) ^ = 0 

An ex 2 unple of zero stiffness can be obtained from a pit extending 
across the tooth profile. 

The equations of motion (21) through (24) contain damping terms 

for all components in the system. For this investigation, the dancing 

in the bearing nearest the driving or driven element has been combined 

with the respective damping of those elements. Damping in the shafts 

[29] 

is due to materiail damping which, by experiment r has been fo\and to 
be between 0.005 and .007. Expressed as a critical damping ratio, a 
representative value of .005 has been assigned for the shafts. For the 
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shafts, the effective dan5>ing is then 


'DS 


DS 



... ( 28 ) 


... (29) 


where 


5 » criticaU. damping ratio 

s 

- .005 


and the shaft masses axe lumped into the masses 

hlso, the effective damping of the gear mesh is: 


CGP^ 



1 

(RBCl)^ ^ 
•^Gl 


(RBC2)' 
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where 


...(30) 


Z- « critical damping ratio of geau: mesh 
G 

System response measurements of geared systems indicated 

that C_ ranges between 0.03 and 0.10. In equation (30) the average 
G 

gear mesh stiffness is used and the equivedent masses of gears 1 and 
2 are concentrated at the base circles to reflect their effect «d.ong 
the line of action. 

The eqpiations of motion (21) through (24) eure numerically 

(281 

integrated using a fourth order Runge-Kutta sdieme.^ This method 
is described in a nxinber of references and will not be repeated here. 

Before the integration can be performed. Initial values, integration 
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time increment cmd integration duration must be determined. For a 
given design condition the initial displacements 6«(0) / 0.(0), 0^(0) 
and 0. (0) are determined by applying the input and output torqtjes to 

Li 

the system. For convenience, the first gear is used as the null point. 
The subsequent driver movement is considered plus and the driven move- 
ment as negative. The initial velocities 0_{O) , 0,(0), 0- (0) and 0_ (0) 

D A L 

have been assigned the anticipated steady state velocities. 

The integration time step must be selected short enou^ to avoid 
inaccuracies and instability in the integration process auad yet long 
enough to minimize computer time. A xneasxire of the optimum time step 
caui be obtained by determining the undaimped torsion^□. natural fre- 
quencies of the system. The undamped eqpiations of motion rewritten 
in matrix form appear below 

[J]{0} + [K]{0} * 0 ...(31) 

\diere the inertia matrix is 



and the stiffness matrix is 



-Ss 

■S>s - '“'ave » 

-■^AVG « »BC1 X BBC2 

0 


0 

-KGP^^^ X RBCl x RBC2 

Ss * “'avg « 

...( 33 ) 
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In equation (33) the weighted average of gear mesh stiffness^ 
KGP^^/ is introduced to simplify the solution for eigenvalues. 
is determined by summing up the discrete stiffness V 2 ilues 7 KGP^, over 
one cycle and dividing by the number of discrete mesh positions in the 
cycle . 

The undcunped eq^oations of motion are solved for the eigenvalues 
and eigenvectors by a Jacobi iteration technique. For the integration 
time step stable solutions have been obtained consistently by using 
one tenth of the shortest system natural period or less than two per- 
cent of the mesh stiffness period. The duration of the integration 
time step is predicated on the time needed for the stairt-up tremsients 
to decay. Review of the output data revealed essentially steady state 
behavior for integration time lengths equal to five times the longest 
system natural period. 

As a first step, the dyneunic force in the mesh as defined by 
equations (25) through (27) is calculated in the dynamic analysis 
subroutine FAST. Next, FAST interacts with the static subroutine 
SLOWM to determine the sharing of the dynamic load, the variation of 
the load through the mesh cycle, the sliding velocity, the maximum 
Hertz pressure and the velocity- Hertz pressure product along the tooth 
profiles. 

The sliding velocity vector relationship at a given mesh position 
can be seen from Figure 21. In vector notation 

SV(k)^ - - RCPl*(k) 0^ - RCP2*(k) ...(34) 

udiere 

RCPl'(k)^ and RCP2'(k)^ are the instantaneous radii to the 
contact point of tooth k at mesh position i. In scalar form equation 


69 


(34) can be expressed as 


SV(k)^ 


/- ‘2 2 


(V^) + (V^) - 2V^V2 cos 


a^) ...(35) 


For the same position the dynamic load QD(k) was estcUslished as 


QD(k). 


KP(k) . 


...(36) 


cuid two dynamic load factors as 

QOT. 

(DFl), * -r— 


...(37) 


(DF2), 


QD(k), 


.•(38) 


where DFl is defined as the ratio of total dynaunic mesh force to the 
total static mesh load. DF2 is the dynamic load ratio for emy given 
pair in engagement. 

Two stress conditions aure evaluated for the contacting teeth: 

1. Hertzian Contact Stress using the equivalent cylinder 
approach for Hertzian deflection outlined in Section 
3.4.6 : 


:^^^i “ y irFA ^RCCPl' 


(k) RCCP2 ' (k) 


H 

where 

RCCPl'(k), BCCP2'(k) » equivalent instantaneous radii 
of curvature 

F • minimum gear tooth face width 

(1 - uj) (1 - uh 
• * . + ^ 


(39) 


E, 


E. 


1 “2 

Using equations (35) emd (39), the product P„(k). x SV(k). 

n X X 
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is determined as an indication of the severity of the 
wear condition at the tooth profile surfaces. 

2. Bending stress of the teeth using a modified Heywood 
formula suggested by Cornell: 


a 


B 


QD(3) . 


X. ,7 


■^COS 0 tl + .26 (^) ] 

j 



TAN 0, 


3S 3 s js 3s 


where 


j - 1 or 2 


For the modified Heywood formula , the position of the maucimum 
stress in the fillet, is fomd by iteration 

- (j+ .016 TAN ...(41) 

where 

t = iteration nximber starting with 1 
2 X 

V ■ -Sff-* 

Y. 

®-4 P * Bn SIN Y • O 
:<c RF^ '3SC 

and the remaining nomenclatTire is as in Figures 23 amd 24. 

Equation (40) predicts the maximum tensile fillet stress within 
about 5% to 10% of finite element methods, It also predicts 
fairly well the location of the peak stress in the fillet. Because 
of its relative ease of use and the expected low stresses in the rigid 
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Figxin 23 - Nomenclat\ire for Modified Heywood Formiila - 

Extemed Geatr Tooth 
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hub and ring, equation (40) is considered to be representative of the 
peak stresses of gear 1 and 2. 
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3,6 COMPUTER PROGRAM 

The development of a digital coirputer program for the comprehensive 
auialysis of the static cmd dynaunic behavior of internal spur gear (ISG) 
drives was one of the tasks of this dissertation effort. The existing 
external spur gear prograun developed at CSU was used as the nominal 
starting point for the structure amd nomenclature of the new ISG drive 
prograun. The ISG computer program package in its entirety, along with 
S2unple output print-out, is included in Appendix C. Highlights of the 
computer program, its structure, nomenclature cmd the input data re- 
quired for the amalysis of an external -internal spur gear drive are 
discussed in this section. 

3.6.1 Program Structure 

The ISG drive con^uter program package is written in Fortran 
IV G1 for use on the Cleveland State University IBM 370/158 digital 
computer. It has been prepared in three modules for operator conven- 
ience (see Figure 25) . Module 1 represents the static analysis dis- 
cussed in Section 3.4. In this module, the operator has the following 
options : 

1. INTERNAL STATIC permits the operator to make changes to 
the prograun and conduct a static analysis. This ability 
to change the program is desirable for future in^jrove- 
roents. The program is compiled after every submittal 
and thus takes over 2 minutes of computer time on the 
CSU IBM 370/158 digital computer. 

2. The cturrent INTERNAL STATIC can be compiled and it 
becomes now EXECUTE INTRSTAT. 'However, before the 
conpiling is done, it is first necessary to delete the 
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Figure 25 - Block Diagram of Computer Package for External -Internal Spur Gears 
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old EXECUTE INTRSTAT. Then the empty space on the file 
is compacted, the old INTERNAL STATIC is deleted, the 
new INTERNAL STATIC version is t reins f erred to the de- 
leted section and then compiled. Now the new EXECUTE 
INTRSTAT is ready. 

3. EXECUTE INTRSTAT permits no changes to the prograun, but 
the computer time is red\K:ed to less than 2 minutes. 

Thus, parametric studies are made most efficiently di- 
rectly from the current EXECUTE INTRSTAT version. 

Option 1 and 3 provide the additional feature to the operator of 
being able to write output data on two tapes. Tape 8 and 9 axe needed 
for the dynamic and stress analysis respectively. 

Nodule 2 represents the dynamic analysis as discussed in Section 
3.5. In this module the operator has the same options as before but 
the uncompiled and conpiled versions are now called INTERNAL DYNAMIC 
and EXECUTE INTRDYNM. Once again, EXECUTE INTRDYNM can be used di- 
rectly for parametric studies. Both versions of the dynamic analysis, 
in addition to input data, require the data from Tape 8, and can write 
dynamic data on Tape 9 for use in the stress cuialysis. 

Module 3 calculates the maximum tooth bending stress at the fillet 
using the stress formulae discussed in Section 3.5. Tape 9 information 
is needed for this analysis in terms of the 50 contact points , and 
whether static or dynauaic loading is used. 

The three modules consist of an executive program and a nusber of 
subroutines. Portions of the executive program are oonmon to all 
modules and some subroutine sections are shared by Modtile 1 and 3. The 
input data is entered into the executive program by use of namelist 
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parameters, and Tapes 8 and 9. Again, some of the namelist parameters 
are common between the three modules* Thus, the resultant overall pro- 
gram package is larger than woxiLd be obtained with a coirbined program. 
However, the advantage of operator flexibility and quick computer tiim- 
around made the "three module" approach far more convenient. Also, an 
initial step has been made towards preparing the analysis package for 
use in the CSU mini-coit^uter HP 1000. 

3.6.2 Description of the Executive Programs and Subroutines 

The executive program of all modxiles is called MAIN. It 
contains the read statements for the namelist parameters and tapes, the 
call statements for the sxabroutines and the write statements for input 
data, diagnostics and output data. Each MAIN program has its own set 
of ncunelist parsuneters, subroutine call statements, and read and write 
statements suitable for the type of analysis of the given module. 

3. 6. 2.1 Module 1 

In addition to MAIN, Module 1 uses MAINl for additional 
executive control and write instructions. The MOD subroutine defines 
the digitized tooth profiles, checks for interference between mating 
teeth and determines the theoretical contact ratio. SLOWM first lo- 
cates the geometric arrangement of the gears with respect to each other. 
Next it establishes the points of contact, number of contacting gear 
tooth pairs, sliding velocity vectors, length of dynamic cycle, and 
loaded contact ratio. Subroutine DEFLECT calculates the deflections 
except the Hertziaui deformation due to a unit load. This unit load 
deflection is combined with the Hertzian deformation in SLOWM to 
determine the stiffness of the individual pairs, the variable mesh 
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stiffness, the static load and the deflection due to this loading, 

3. 6. 2. 2 Modiile 2 

The FAST subroutine of this modvile analyzes the three degree 
of freedom, four mass, mathematical model of the geared torsional sys- 
tem depicted in Figure 22. Calculations in FAST are based on a dynamic 
cycle which starts with the initiation of contact on a tooth entering 
the contact zone and ends with the initiation of contact with the tooth 
following it. The length of this cycle is established in subroutine 
SI/DWM by examining the stiffness fvinctiai. The position of tooth no. 3 
when tooth no. 4 comes into contact* is defined as lEP. Consequently, 
(IEP-1) is the endpoint of the dynamic cycle started when tooth no. 3 
came into contact (Figure 27) . FAST calculates the dynamic force in 
the mesh as defined by equations (25) through (27) . The integration is 
performed in two small subroutines using the very efficient Runge-Kutta 
integration scheme by Freuiks. ^ ^ The RKUTTA subroutine contains the 

integration step size and keeps track of the iterations across the 
integration interval. The actual integration is done in subroutine 
MORERK. Subroutine VIES uses a Jacobi iteration technique to determine 
eigenvalues of the gear train. This information is returned to FAST 
for determination of the natural frequencies, eigenvectors, integration 
time step and duration. FAST also calculates the instantaneous angular 
position and velocities, sliding velocities, Hertzian pressure, dynaunic 
loads, dynaunic load factors, and transmission ratio. Subroutines STORE 
auid XPLOT are used to store the data and then plot as many as four 
dependent variables against a single independent variable (see Appendix 
C for sample- output plots) . 
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3. 6. 2. 3 Module 3 

Module 3 is set aside for calculating the maximum bending 
stress of the tooth fillet. Portions of MAIN and MAINl have been 
modified for executive control and printing of pertinent information. 
Subroutine MOD has been modified to provide not only the tooth profile 
points but also a more refined breakdown of the fillet contour. This 
additional refinement is needed for the iterative seaurch routine of 
the maximum fillet stress location as indicated by equation (41). The 
actual stress calculation is done in subroutine CORNEL. 
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CHAPTER IV 

RESULTS, DISCUSSION AND SUMMARY 

4.1 RESULTS AND DISCUSSION 

4.1.1 Introduction 

The internal spur gear (ISG) analysis methodology, which was 

developed for this report, was used to perform a series of con^arative 

parametric studies in order to assess the ISG drive performance. Since 

there are no ISG drive performance data available, comparisons were 

made with known solutions of external spur gear drives. In particular, 

[191 

the results of the ESG study by Kasuba and Evans can be used to 

compare the static and dynamic performance of internal versus external 

spur gear drives under identical load, speed and geometry conditions. 

From the study by Kasuba and Evams, this author selected one external 

gear set ^ich was of practical interest to the internal qeax design. 

The selected external spur gear (ESG) set consists of the following 

design paLrameters: 

Nuinber of Teeth on 
1st and 2nd Gear 

Pressure Angle 

Diametral Pitch 

Face Width 

The analyses of the ESG set consider variations in tooth profiles, 
tooth support stiffness, critical damping ratios and shaft stif&iess* 
For the ISG drive, static and dynamic computer analyses were conducted 


32 & 96 

14.5® 

8 

25.4 mm 
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under identical conditions. Next, ISG drives of practical interest 
were investigated for the same output information to allow comparisons 
between external and internal drives. Additional analyses were con- 
ducted to investigate the effects of radial deflection of the bearings, 
shafts, hub and gear rings on the ISG drive performance. Also, the 
tooth bending stress program was used to determine the static and 
dyneunic stress for all fifty gear mesh positions. 

4.1.2 Static Analysis 

The dynaunical model shown in Figure 22 was also used in the 
static analysis (with all masses taken to be zero) , The behavior of 
the internal and external geau: teeth were investigated for various 
structTire support conditions and load magnitudes. In the static analy- 
sis the information of interest for the gear teeth consists of the 
deflection, stiffness, load sharing, bending stress, unit sliding 
velocities and Hertz pressxire. 

4. 1.2.1 Comparison of ISG and ESG Set Performance 

The results presented in Table 1 and Figure 26 show a com- 
parison of the mesh stiffness chauracteristics between error-less 
external and internal spur gear sets. The results indicate the in- 
fluence of the tooth support stiffness on the overall gear mesh stiff- 
ness for both ESG amd ISG sets. By increaising the support torsional 
stiffness (higher HSF, Appendix B-3.3), the loaded contact ratio 
decreases, mesh stiffness increaises, changes in instauitameous trans- 
mission ratio decrease, and sensitivity to gear tooth errors increases. 
The opposite occurs by decreasing the support stiffness. In the case 
of the ISG set, the geau: ring stiffness can be reduced effectively only 
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by decreasing the thickness of the rim (see Figures 3 and 4). 

Table 2 demonstrates the significantly higher theoretical contact 
ratio of the ISG versus the ESG set. Also, the tcibulated results in 
Table 1 indicate siibstantial changes in contact ratio with increasing 
loads and/or support flexibility. For different load magnitudes ranging 
from 88 to 700 N/m (500 to 4000 lbs. /in.) the contact ratios change by 
13% and 7% for the respective ISG and ESG sets. At the maximum load con- 
dition, the contact ratio of the ISG set is 3.09, i.e., three pairs are 
in continuous contact. This high contact ratio is obtained by the 
internal gear set, even with high gear tooth support stiffiiess (HSP-1.0). 

Profile errors and pitting affect the mesh stiffiiess characteristics 
of the ISG and ESG sets to varying degrees. For the case of only the 
pinion or gear having a narrow pit 0.5 mm wide (0.02 in.) at the pitch 
line, the torsionally stiff ISG gears absorbed the fault without a 
change in the mesh stiffhess pattern whereas the ESG gears could not 
absorb this fault. Thus, causing significant interruptions of the mesh 
stiffness pattern. Only v4aen the pit width was increased to 2,0 mm 
did the error affect the ISG mesh stif&ess characteristics (Figure 27) . 
The unabsorbed error caused non -contacting zones with resulting sub- 
stcuitial changes in the mesh stiffness characteristics, i.e., the 
flexibilities in the mesh were aible to bridge the non -contact zones. 

For the pit shown in Figure 27, the flexibility of the mesh was able 
to absorb a portion of the error by eliminating about 60% of the mesh 
stiffness interruption. The reason for this effective bridging of 
the non-contact zones can be found in the torsional ring flexibility 
and the inherent high contact ratio of the ISG set. In normal contact 
ratio gears as can be encoxmtered with ESG sets, the mesh stiffcess. 
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TABLE 1 

EFFECTS OF GEAR HUB/RING FLEXIBILITY ON MESH STIFFNESS, 
TRANSMISSION RATIO AND CONTACT RATIO FOR 
INTERNAL AND EXTERNAL SPUR GEAR SETS 


Gears; 32 & 96 T, 8 DP, 14.5“ PA, CR.J, = 2.625 (INTERNAL), 
CR^ =2.14 (EXTERNAL) 



Normal 

Load: 

4450 N (1000 lb. 

) or 175 N/mm 

(1000 

lb. /in.) 


GEAR 

TYPE 

RHlf 

mm 

mif 

mm 

^^max 

N/m 

KGp 

N/m? 

HSF 

ATR* 

% 

CR 

ISG 

12.7 

218.5* 

6.12 X 10® 

2.41 xlO^ 

.88 

1.8 

2.853 


12.7 

218.5 

6.15 X 10® 

2.42 X 10^ 

.89 

1.8 

2.853 


12.7 

171.7 

6.15x10® 

2.42 xio"* 

.89 

1.8 

2.853 


47.2 

156 

6 .94 X 10® 

2.73x10^ 

1.0 

1.5 

2.81 

ESG 

10 

14.5 

3.07x10® 

1.21 xlO^ 

.476 

2.4 

2.47 


12.7 

18.3 

3.80 X 10® 

1.50 xio'* 

.591 

1.9 

2.42 


12.7 

38.1 

5.08x10® 

2.0 xio'* 

.794 

1.6 

2.36 


38.1 

114.3 

6 . 36 X 10® 

2.5 xlO^ 

.992 

1.0 

2.32 


47.2 

148.8 

6.45x10® 

2.54 X 10^ 

1.0 

1.0 

2,32 


KG = maximum attainable stiffness in meshing arc 
max ^ 

= FHl = FH2 = 25.4 mm (1.0 in.) 

RRCl = 47.25 mm RRC2 = 156 mm for internal gears 

RRCl = 47.25 mm RRC2 = 148.84 mm for external gears 


KG 


KG^ 


max 


All gears without errors or modifications 
*Rim thickness of ring gear = 10 mm 
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TABLE 2 

LOAD EFFECTS ON MESH STIFFNESS, TRANSMISSION RATIO 
AND CONTACT RATIO FOR INTERNAL AND EXTERNAL 
SPUR GEAR SETS 


Gears: 32 & 96 T, 8 DP, 14.5“ PA, CR.p = 2.625 (INTERNAL), 

CR^ = 2.14 (EXTERNAL), HSF » 1.0 


GEAR 

TYPE 

LOAD 

N/m 

KGmax 

N/mm 

KGp 

N/mm2 

atr 

% 

CR 

ISG 

88 

6.94 X 10® 

2.73x10^ 

1.15 

2.731 


175 

6.94 X 10® 

2.73x10^ 

1.51 

2.81 


350 

6.94 X 10® 

2.73 xlo"* 

1.83 

2.895 


525 

6.94x10® 

2.73 X 10^ 

2.05 

3.02 


700 

6.94 X 10® 

2.73 xlO^ 

2el6 

3.09 

ESG 

88 

6 . 36 X 10® 

2.5 xlO^ 

0.8 

2.29 


175 

6.36x10® 

2.5 xlO^ 

1.0 

2.32 


350 

6.36 xlO® 

2.5 xlO^ 

1.0 

2.38 


525 

6.36 xlO® 

2.5 xlO^ 

1.8 

2.43 


700 

6 . 36 X 10® 

2.5 xio'^ 

2.2 

2.45 


KG “ ma ximum attaixiable stiffness In meshing arc 
max 

■ FHl ■ PH2 ■ 25.4 mn (1.0 in.) 

RBCl m 47.25 nsD RRC2 " 156 nm for intemad gears 

RRCl “ 47.25 mm RRC2 • 148.8 smi for external gears 

KG 


KG_ 


max 

F 


7 


All gears without errors or modification 




Gears as listed in Table 1; HSF = 1.0 and Load = 175 N/m 
Angle of approach is negative; Angle of recess is positive 









7.0 I r WITHOUT PIT 


86 



8 - 


01 X (ui/N) oyL *sszmdijs Hsaw 
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01 X (m/N) HS3N 


Figure 28 - Effect of Sinusoidal Error on ISG Mesh Stiffness 
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KGr wotild equal zero when these non-contacting zones due to errors occur. 

Sinusoidal errors of 0.013 mm (.0005 in.) on either the pinion or gear 

of the ISG drive caused non-contact zones and the resulting changes in 

mesh stiffness characteristics of Figure 28 due to its high torsional 

rigidity. In contrast, a torsionally flexible ESG drive was able to 

. . [19] 

cdjsorb the fault without affecting the mesh stiffness characteristxcs. 

The analysis procedure can also be used to investigate other sur- 
face fault-error combinations acting on both gears. For example, errors 
of Figure 16c caused mesh stiffness reductions over a longer span of the 
meshing cycle than were obtained when the identicaO. error was on one gear 
only as in Figure 28. Thus, it must be concluded that each profile con- 
dition and mesh system flexibility will cause unique mesh stiffiiess 
patterns. The geau: tooth contacts due to deflections do not coincide 
with the theoretical line of action. This resxilts in noninvolute action 
producing variations in the transmission ratio, ATR. For the investi- 
gated cases in Tables 1 amd 2, the trend in ATR for the ISG amd ESG sets 
was the same. The magnitude of the ISG ATR was about 50% higher but not 
exceeding 2.7%. However, the transmission ratio of the ISG set is higher. 

4. 1.2. 2 ISG Drives of Practical Interest 

Review of the maximum stiffness, KG , in Table 1 reveals 

max 

a 13% higher stiffness value for the torsionally stiff ISG drive which 
does not rediice significantly with changes in hub or ring radii changes. 
The reason for this small vauriation in stiffness can be attributed to 
the rather rigid internal gear ring. In order to reduce this circian- 
ferential stiffness of the ring, the rim thickness has to be reduced to 
20% or less of the face width of the gear teeth. However, ISG drives 
of practical, interest do not exhibit such smaill rim thicknesses in order 
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to minimize deflections in the radial direction (Figures 3 and 4) . The 
inference is that ISG drives will be inherently stiffer than ESG drives 
and, because of their closely matched contours, will be less tolerant 
to sinusoidal imperfections or similar irregularities of the surface. 
This intolerance to surface irregulcirities can become limiting since 
grinding of the internal teeth is not practical. 

4. 1.2. 3 Effects of Radial Deflection on Static Performance 

Radial deflections due to bearings, shafts, pinion hub and 
gear ring cause radial movement of the gears in the rotating plane with 
a result 2 uit reduction of contact length of the gears. This fact CeUi 
readily be deduced from a plot of contact ratio versus radial deflec- 
tions (Figmre 29) . The plot indicates a sharp reduction of contact 
ratio with small radial deflections. At .025 mm (.001 in.) radial 
deflection, the reduction of contact ratio levels off and reaches the 
theoretical contact ratio, CR^^, of 2.625 at a radial deflection of 
0.05 mm (.002 in.). AdditioneJ. radial deflection causes the contact 
ratio to drop below the CR^. For the gear auid loading combination of 
Figure 29 practical combined bearing, shaft, pinion hub and gear ring 
deflections can be held to 0.05 mm (.002 in.). Thus, it appecurs that 
the actual loaded contact ratio probably will not exceed the CRj^, because 
of the ever present radial deflection unless design and assembly prac- 
tice offsets the radial deflection with an equal amount of reduction in 
center distance. It is interesting to note that the combined circiim- 
ferential deflection of the two ISG drive gears of Figure 29 has a 
significemtly greater auid opposite effect on contact ratio than does 


the radial deflection. 




CXJMPACT RATIO 
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H2VDIAL DEFLECTION, mm 


Figure 29 - Effect of Radial Deflection of 
ISG Drive on Contact Ratio 
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For the load condition of Figure 29, the two deflections offset 
each other at a radial deflection of .05 mm (,002 in.) and a corres- 
ponding maximum circumferential deflection of .015 mm (.0006 in.). 

The contact search for radial deflections of 0.025 mm or less 
can be accomplished with the contact search method as discussed in 
Section 3.4.4. For larger radial deflections the contact search 
method was changed to include non-uniform rotation ("wiggling") for 
establishing of the limiting mesh points. 

4.1.3 Dynamic Analysis 

The geeu: train of Figure 22 was modeled for the dynamic analysis. 
The solution of the dynamic equations of motion (equations 21 through 
24) leads to dynamic loads which depend on the magnitudes of the mass 
moments of inertia of all elements, shaft stiffness, transmitted loads, 
gear mesh stiffness characteristics, damping in the system, amount of 
backlash and speed. 

4.1. 3.1 Comparison of ISG and ESG Drive Dynamic Performance 

For comparison purposes, identical geometry and operating 

conditions were applied to the ISG drive model of Figure 22 and the 
[19] 

ESG drive model. Figures 30 and 31 show the results of a series 

of analyses on error-less ISG and ESG gear trains using low and high 
stiffnesses for the shafts cuid teeth supports and varying amounts of 
damping. Review of the resulting ciirves indicates a smoother per- 
formance in terms of lower peak dynamic factors for the ISG drive. 

Figtire 32 demonstrates the dramatic effect of severe mesh stiff- 
ness interruptions due to sinusoidal errors on dynamic performance. 

For this type of error the two drives exhibit similarly violent 
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dyn£unic fluctuations. The performance of the ISG versus ESG drive due 
to the effect of a surface pit is considerably smoother (Figure 33). 

In this error condition, locating the pit on either the pinion or gear 
of the ISG drive had little effect on performance. 

There is a requirement for a minimum amount of damping to prevent 

rgi 

the Mathieu-Hill type instabilities. For the considered ISG drive 

of 32 and 96 gear teeth, there was no such instability encountered 
within the operating ramge of 0 to 12,000 rpm with “ *05 and Cg = .005. 
Because of the large support bearing of the ring gear (Figures 3 and 4) 
the conbined bearing and gear mesh damping, (q, is probably neeu: .15, 
and thus Mathieu-Hlll instabilities are unlikely to be encountered with 
ISG drives. 

The XS6 analysis procedure has the eaped^llity for analyzing the 

distribution of the dynamic loads, dynamic factors (Figures 30 through 

33 ) , load sharing, contact Hertz stress (P ) , contact stress-sliding 

n 

velocity product (FV) and maximum tooth bending stress for the entire 
meshing zone. Figures 34 and 35 show the range of the maximum and 
maximum PV values corresponding to the dynamic conditions illustrated 
in Figures 30 and 33. In general, the higher contact ratio ISG drive 
appeared to show lower values. 

Pigture 36 shows the maocimim tensile tooth bending stress versus 
all fifty gear mesh positions for the static and the 8000 rpm operating 
condition of the ISG drive depicted in Figure 30, case d. As expected, 
the pinion shows higher stress values both for the static and dynamic 
condition. The peak stress values occur in the gear mesh position 
range from 14 to 18 and 33 to 38. These geau: mesh positions represent 
the change over from three to two tooth contact. Figures 37 and 38 show 
the stress pattern for the ESG drive. 


93 


o 

o 


- o 

O 'O T3 
O tTj (N m 



m CM (N iH iH O 


CT» 


00 




VO 


in 




CM 


O 


§ 

cz: 


u 

s 

u 




0e« 

Q 

Q 

£e 3 

a< 

CO 


'H0J.0VJ oiwvj^a 


Figure 30 - System Flexibility Effects on Dynamic Factors 
of Error-Less ISG and ESG Drives 
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Figure 32 Influence of Sinusoidal Error on Dynamic Factors for 

ISG and ESG Drives 
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FigurG 36 - Maximum Tooth Banding Stress Versus Gear Mesh Position 
for the Static and 8000 rpm Operating Condition of the ISG Drive 
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Figure 37 - Maximum Tooth Bending Stress Versus Gear Mesh Position 
for the Static and 8000 rpm Operating Condition of the ESG Drive 
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Figure 38 - Maxiimim Tooth Bending Stress Using AGMA Formula Versus Geao: 
Mesh Position for the Static Operating Condition of the ESG Drive 
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4.1. 3.2 ISG Drives of Practical Interest 

For similar construction, ISG drives are torsionally stiffer 
than ESG drives because of the large outside diameter of the internal 
gear. Thus, the torsionally soft case of HSF = .65 and mass moment 
of xnertia, = 0.017 m -kg (Figures 30 through 33) , can only be 
achieved with a severe reduction of the internal geeu: rim thickness 
(Figures 3 and 4) . The dynamic load factors of ISG drives of prac- 
tical interest with various mass moments of inertia, J , and differ- 
ent gear mesh and shaft stiffnesses are shown in Figure 39. For the 
four cases investigated, the best performance, in terms of lowest 
gear dynamic load factors, was obtained for the ISG drive of high 
mass moment of inertia, high shaft stifftiess and low torsional stiff- 
ness, HSP ■ .88. The ESG drive of equal shaft stifftiess, HSF - 1.0 
(Figure 30, Curve c) exhibited the highest dynamic load factors. The 
results of the analysis of practical ISG drives demonstrates the need 
for tuning of the various constituent elements in the model of Figure 
22. The general trend indicates smoother performance from ISG drives 
of practical interest versus similar ESG drives. 

4. 1.3. 3 Effect of Radial Deflection on Dyneunic Performance 

Figure 40 shows three cases of ISG drives subject to different 
radial deflections. Review of the results shows larger dynamic load 
factors with increase in radial deflection. The performeuice of the 
ISG drives exhibits considerable missing and "badehitting" of the 
gear teeth similar to the performance of geeurs with sinusoidal profile 
errors. Significant reduction of the dynamic load factors was evident 
with increase in gear mesh damping. 
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Figure 39 - Dynamic Load Factors for Various ISG Configurations 
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4.2 SUMMARY AND CONCLUSIONS 

A new methodology has been developed for the static and dynamic 
load and stress analysis of the internal spur gear (ISG) drive. Prior 
to this report, there were no established methods for the above analy- 
ses. The currently published design techniques for ISG drives reflect 
the technology of the 1950' s. 

The analysis procedure is applicable to involute profiles cind 
minor deviations from this profile as a result of modifications/ imper- 
fections and circumferential deflections. Because of the potential 
noninvolute profile and the effect of radial deflection on tooth posi- 
tion, an iterative procedure is used to calculate the statically 
indeterminate problem of multi-tooth contacts, circumferential deflec- 
tion and contact ratio. The developed method can be used in gear 
combinations leading up to and exceeding the "very high contact ratio" 
(VHCR) of three. The static anad.ysis can also be adapted for deter- 
mining the gear mesh stiffness of a planet and ring gear ass^obly. 

For this adaptation it is necessary to supply the ring gear stiffness 
relying on strength of materials, finite element or experimental mecms. 
The maximum tooth bending stress of the external and internal gear is 
determined by "Cornell's method". This method is a modification of 
the empirical foirmula for stress of loaded projections based on photo- 
elastic experiments by Heywood. 

The new methodology was compiiterized with the computer package 
consisting of three modules which perform the static, dynamic and 
stress analyses respectively. The output from the three modules in- 
cludes the static and dynamic loads, variations in transmission ratio, 
sliding velocities, maximxim contact pressures and tooth bending stress 
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acting on the gear teeth at the fifty mesh positions. The resiilts of 
the parametric con^ter studies yielded the following conclusions: 

1. For equal geometries , the ISG drive achieves higher 
contact and transmission ratios than the ESG drive. 

For the investigated 32 and 9^ tooth combination, the 
heavily loaded ISG drive reached the "very high contact 
ratio" of three. The unloaded or theoretical contact 
ratio is 2.625. 

2. Because of the high contact ratio and the concave- 
convex tooth contours, the ISG drive is able to absorb 
pits across the surface better than ESG drives* How- 
ever, the opposite effect is encountered when sinusoidal 
or similar surface irregularities are present. 

3. ISG drives tend to be torsionally stiffer and heavier 
than ESG drives because of the large intemaU. gear. 
Nevertheless, when both systems are sized optimally, 
then the ISG drive performance in terms of peak dynamic 
loading is better. 

4. Radial deflections of shafts, bearings and gears reduce 
the contact ratio of ISG and presiamably ESG drives. 

The effect on mesh stiffness is minor but the pattern 
is chamged. The dynamic loading increases with radial 
deflection. 

5. The dynaunic factors of ISG drives cam be reduced with 
increases in damping ratio. This effect is also evident 
in drives experiencing radiaU. deflections. 

The maximum product of Hertz stress aind sliding velocity 


6 . 
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is consistently lower than the ESG drive. The naximum 
Hertz stress performance is distorted because of the 
different peeUc dynamic loading of the various cases that 
were investigated. In general, the maximum Hertz stress 
appears lower also for the ISG drive. Peak bending stress 
of the pinions were the same whereas the gear of the ISG 
drive experienced an 18% lower stress value. 

The derived new euialysis procedure has established a method 
exclusively for the emalysis and prediction of dynamic performauice of 
the intemal-extemad. spur gear set. The list of advantages, as de- 
rived from this study, of the ISG drive over the ESG drive is quite 
isqpressive, and should lead to renewed interest in applying the ISG 
drives in advanced transmissions. 

Recommendations for fvrtiire work might include a finite element 
method for the ring geeu: deflections and an efficient contact search 
technique which would satisfy the intentioncLLly built-in large radial 
deflections such as foisid in the planetary gear rings. 
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APPENDIX A 

SPUR GEAR FORMULAE AND INVOLUTE PROFILE DEVELOPMENT 

This Appendix lists st£uideu:d spur gear geometry relations and 
develops the involute profile for the internal and external tooth. 
It forms the basis for the definition of the actual tooth profiles 
in Section 3.4.3. 

A-1 Stamdaurd Spur Geau: Relations for the ISG Drive 
A-2 Development of the Involute Profile 
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A- 

1 STANDARD SPUR GEAR 

RELATIONS ^ FOR ' 

THE ISG DRIVE 



Definitions 


P 

- diametral pitch 

F 

= face width 

M 

= module 

B 

= backlash 

r 

* pitch radivis 

*n 

= normal pressure angl 

r 

o 

= outside radius 


= pressure angle 

r 

r 

= root radius 

0 

= involxite polar angle 

r 

a 

= addendum radius 

Y 

= angle between pitch ; 




and center of tooth 


* dedendum radius 





a 

= addendum 


= base circle radius 





d 

= dedendiim 


= limit radius 






= whole depth 


= fillet radius 

u 




e 

* roll angle 

""t 

= edge raditis of 




generating tool 

u 

= interval of contact 

N 

“ number of teeth 

“g 

* gear ratio 

C 

= center distance 


* contact ratio 

P 

= circular pitch 

» ir^ut torque 

Pb 

* base pitch 


* tcuigential load 



Subscripts 


1 

=* external gear 

o 

= outside 

2 

“ internal gear 

r 

* root 

b 

■ base 

t 

= tangent, total 

F 

- fillet 

T 

» tool 

G 

« gear 

P 

» pitch 

1 

* limit 

in 

= input 

n 

» normal 
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Standard Fonn\ilae 


2P 


^2 " ^1 


Pb * P ♦n 


(h^ + a - r^)' 


Rj = 0.7 (r^ + j, ^ ^ j. j 


r - r COS 4^ 


--1 ’^b 


♦ - COS — 

a 

- -(r^ - a^) SIN ^ SIN ^ 

"2 " *1 " ^n 


Tj SIN 4,^ > 


'2 ^ '^2 


* COS 


\ “ cos 


\ ^2 


e * TAN ()) 


A-1 

A-2 

A- 3 
A-4 

A-5 

A-6 

A-7 

A-8 

A-9 

A-10 

A-11 

A-12 

A-1 3 

A-14 

A-15 
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0 = TAN (|) - 

n n 


T. 

xn 



hfc = a + d 

r = r + a 

o 


r 

r 


r - 


d 


ro - r - a 
« r + d 

M - 2Ld. 
M - — 


A-16 

A-17 


ECTERNAL 

GEAR 

A-18 

A-19 

A-20 

EXTERNAL 

GEAR 

A-21 

INTERNAL 

GEAR 

A-22 

INTERNAL 

GEAR 

A-23 



A-24 
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DEVFTfQPMENT OF THE INVOLUTE PROFILE 

The construction of the involute spur gear tooth follows exact 
geometric relations as indicated by the standard formulae in Section 
A-1. Thus, from knowledge of a few parameters all the other parameters 
can be determined. General design practice starts with an assumption 
of the number of teeth, diametral pitch, addendian, working depth, fillet 
radius and pressure angle at the pitch point. The radial distance, R, 
to any point along the involute line is now determined from the base 
radius and the pressure angle at that point (see Figure A-1) . The in- 
volute can then be drawn by connecting finely spaced radial points with 
straight lines. In the computer program either one or two hundred 
points are used depending on the size of the tooth. Any position be- 
tween two points is determined by linear interpolation. The external 
and internal tooth have the same involute profile. The distinction 
between the two arises from the fact that their tooth center is either 
on the concave or convex side of the involute profile. Fig\ires A-2 and 
A- 3 respectively show the construction of the internal and external 
involute tooth along with the fillet configuration and the position of 
the local and global tooth coordinate system. 
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Figure A-1 - Involute Development 
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itemal Gear Tooth Involute Profile 
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APPENDIX B 
DEFLECTIONS 


The radial and circumferential deflections of the external -internal 
gear support bearings, gear ring and teeth as discussed in Section 3.4.6 
are summarized in this Appendix. 

B-1 Bearing Deflections 

B-2 Radial Ring Gear Deflections 

B-3 Circumferential Deformation of Gear Teeth 

B-3.1 Deflection of Point of Contact Due to Deformation 
of Teeth 

B-3. 2 Deformation of the Teeth Due to Rotation of Their 
Foundation 

B-3. 3 Deflection of the Teeth Due to Circumferential 
Deformation of the Rim and Gear Ring 


B-3. 4 


Hertzi 2 ui Deformation at Contact Point 
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B-1 BEARING DEFLECTIONS 

Rolling element bearings are preferred for their use in ISG drives 
because of their load carrying capacity, durability and low maintenance. 
The deflection in rolling element bearings is mainly due to Hertzian 
contact deformation. Because the maximum elastic contact deformation 
is dependent on the rolling element loads, it is necessary to cuicdyze 
the load distribution occurring within the bearing prior to determina- 
tion of the bearing deflection. 

. [26] 

Harris suggests an approximate solution for the maximum ball 
or roller load, as: 

5 

" Z COS a • • • 

where 


= radial load 

Z * number of balls or rollers 
a = contact auigle 

The deflection of a ball or roller bearing are respectively: 


6 

r 


1.58 X 10 


-5 


2^x 


COS o 


6 

r 


» 4.33 X 10 


33/4 

max 


4 COS a 


. . . (Bl-2) 


... (Bl-3) 


where 


D = ball diameter 
t ” length of roller 
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B-2 RING DEFLECTION 


The ring gear of Figure 3 is subject to a single load which is 
reacted by an inclined roller bearing. The ring gear consists of the 
internal teeth, intermediate tooth support ring, outer ring and cylin- 
drical support. Under load the intermediate and outer ring act rigidly 
and thus transfer the load uniformly into the support cylinder. This 
represents nearly cauitilever loading and its deflection represents a 
"best possible" solution. In the extreme, the cylinder can be looked 
at as reacting to the equal and opposite loading of a ring. This 
deflection represents the "worst possible" condition. 

The applicable equations for either ring deflection are: 

Cantilever 


6 

r 


!riL 

3 El 


. . . (B2-1) 


\^ere 

I = — {U - d ) 

64 

D = outside diameter of cylinder 
d = inside diameter of cylinder 


Rinc 


6 

r 


.0745 


F R 
r 

El 


3 


... (B2-2) 


... (B2-3) 


where 

R » mean diameter of cylinder 
I * moment of inertia of ring cross-section 



12 

h « ring thiclcness 


b =« ring gear width 
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B-3 CIRCUMFERENTIAL DEFORMATION OF GEAR TESTH 

B-3.1 Deflection of Point of Contact Due to Deformation of Teeth 
Weber solved for the deflection at the point of application 
of an external gear tooth by equating the stress energy to the deform- 
ing work, ^ In his formulation, the stress energy is composed of 

the partial energies due to the bending moment, the shearing force and 
the normal forces as seen in Figure B3-1. 


- 06 » - — dv + - — d + i. si d 

2 ® 2 El * 2 GFfi 2 EA 


... (B3-1) 

In Figure B3-1 the applied force Q is transferred to the center of the 
tooth resulting in an equivalent loading set 


M - Q(Yg - Y^) COS 6 ...(B3-2) 

V « Q COS 0 ...(B3-3) 

N = Q SIN 9 ... (B3-4) 

and for this cantilever loading 

I = ...(B3-5) 

A = 2x ... (B3-6) 

F = |- for spur gear teeth . . . (B3-7) 


if we replace, E by 


2(1 + u) 

and y » 0.3 

then siabstitution and sin5>lification leads to 
deflection at the point of contact 


... (B3-8) 


the expression for 


5 = ^ COS 9^ 


Y (Y - Y)^ 

[12 / ° — — dy + 3.12 

° (2X)'’ 


(1 + 


TAN^ 9. i^c d^, 
3.12 ’ o 23? 


. . . (B3-9) 
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In this equation x and y are the local coordinates of the involute 
or modified involute profile. By considering the digitized profile 
points of Figure B3-1, it is possible to solve for the deflection at 
the contact point by numeric integration. In this investigation the 
integration is carried out to the intersection of the fillet with the 
root radius, RRCl. Thus, a slight improvement over Weber's disregard 
of the fillet has been achieved. 

Comparison of the more pronounced profile curvature of the external 
versus the internal gear tooth leads to the conclusion that Weber's 
solution equad.ly applies to the internal tooth (see Figures B3-1 emd 
B3-2) . Equation B3'>9 can be used for the deflection of the external 
or interned gear tooth by substituting the appropriate local coordinates 


for the respective tooth. 




Figixre 
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B-3,2 Deformation of the Teeth Due to Rotation of Their Foundation 
Weber investigated the effect of the elastic support of the 
teeth on the deflection at the contact point. For this investigation 
he assumed that a rigid tooth is acting on a semi-infinite support struc- 
ture represented by the gear hub or ring, and that the loading is trans- 
ferred to the support as indicated in Figures B3-3a through B3-3d. He 
then proceeded to find a stress fxmction which satisfied all of the 
boundary conditions for the semi-infinite support. Again, equating the 
deforming work auid stress energy: 

j Q<5 = + S2 *^33 * •••(^3-10) 

where ^22 ^33 factors whose determination is outlined 

as follows. 

By potential fmctions we obtain the deflection in the Y-direction, 


where 


2(1 


boundary 


-LL 


6M 


irb ^2 


r 2 ,b»2, a |b + X| xb, 

{X - (J) } 1-— I - — ] 


. . . (B3-11) 


and 


b = 2 X . 

mm 


X . 

mm 


... (B3-12) 

x-dimension at fillet-to-root intersection 


Now consider the work done by the load at the bomdeury due to M; 

2 . 


the expression for the strain energy M is 


Ci^M 


X r+b/2 » , 

2 -b/2 ^boundary boijndaury 


... (B3-13) 


vdiere ay is the stress. Therefore, 


<3y, 


botmdary 


2x 6M 


By substituting 

Dounoary boundary 


... (B3-14) 


in equation (B3-13) cuid integrating 
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^ m2 - (1 - U ) , 6M , U , 

irE ^ 2^ * . 2 * ^ 16^ 
ijd d 


. . . (B3-15) 


C _ 9 (1 - u^) 1 9(1 - u^) 

11 TT E ..2 2 


. . . (B3-16) 


By similar procedure we get 


^ (1 + u) (1 - 2y) 
2E£b 


. . . (B3-17) 


^22 ~ 


2.4(1 - u ) 


. . (B3-18) 


c . c (1 + ^ 

^33 22' 3.1 ' 


. « (B3-19) 


Substituting the egressions for ^12' ^22 ^33 ^ equation 

(B3-10) and the load set M - Y Y COS 6, V » Q COS 9 and N - Q SIN 6 

c 

we get: 

The deflection of the point of contact dtie to rotation of the 
support structure 


O 2 ^*2 ^c TAN^ 6 

6 = ^ COS 0 [ ^ + -r^ + 1.4 (1 + , ■■■■ - ) ] ... (B3-20) 

£f , z O 3.1 

b 



Figure B3-3 - Stress Distribution on Support Structure 
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B-3.3 Deflection of the Teeth Due to Circumferential Deformation 
of the Riig and Gear Ring* 

Assume the external gear is held rigidly as in Figure B3-4. Under 
the effect of tcuigential force Q, the line HW is deformed and takes the 
position HW* . Any differential element ABCD is deformed to A*B*C*D*. 

The iftovement of A*B* relative to CD* gives an ang\ilar rotation d6 about 


the center of the gear. 

Total angular displacement of the point W =» 0^ » /d0. From Figure 
B3-4, rd0 * ydr; therefore 

de--— ...(B3-21) 

G r 

and 0^ - / ® r ^ • • • (B3-22) 

u r^ G r 

For eqailibrluBir the total shearing force on any concentric sur- 
face is equal to the applied torque T; therefore 


T ■ 2ir r Ft 

T ,^o to JL_ _1_. 

®t “ 2 itLG •'r. 3 “ 4 itLG ' 2 " 2' 

1 r r. r 

1 o 


. . . (B3-23) 
... (B3-24) 


The deflection, 6 , of the point of contact due to the circumfer- 
E 

ential deformation due to force Q 


6 


E 


Q COS 0 2 1 . 

47tFG ^i ' 2 2 ^ 

r., r , 
il ol 


... (B3-25) 


\^ere 

F » hub face width 

* radius to the contacting point 
r^l * outside hub/ria radius 

- effective radius of circumferential hub fixity 
*These derivations are primarily based on [14]. 
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For the internal gear ring deflection a similar rationale as for 
the external gear leads to the expression for deflection (see Figure 


B3-5) . 


Q COS 9 2 

4ttFG ^2 



.. . (B3-26) 


vrtiere 

r^ = radius to contacting point on internal gear 

r ^2 ® root circle radius of gear ring 

r - = effective radius of circumferenticd gear 
o2 

ring fixity 

A torsionally rigid hub is obtained if the effective h\ib fixity 
radius coincides with the root circle , i.e. , 6 = 0. The hub stiffness 
factor, HSF, is used to indicate a degree of influence of the hub/ring 
flexibility on the overall gear mesh stiffness. 


HSF 


KG 


max 


KG 


s max 


. .. (B3-27) 


where 


KG = maximum mesh stiffness with torsionally 

s max 

rigid hubs/rings; 6=5=0 

£j ± 

KG = maximum mesh stiffness with designated 


max 


hubs or rings, 6 euid 6_ 0 


A combination of a rigid external gear hub and rigid internal 


gear ring is identified as HSF = 1. 
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B-3.4 Hertzian Deformation at Contact Point 

Weber applied Hertz's solution between two loaded cylinders in 
contact to the deformation of spur gears. In his formulation he used 
the following rationale: 

The force at the contact point is distributed to the center of 
the tooth and then transmitted to the gear body. The tooth is can- 
tilever v^ich has an equivalent loading at the center consisting of 
bending r compression and shear. The shear load is distributed over 
the cross-section in the form of a parabola. The Hertzian stress is 
dxJB to the shear ooroponent and reaches to the center of the tooth vdiere 
it is transmitted as a transverse stress. The Hertziam compression 
is calculated from the point of contact in the direction of the applied 
force to the center of the tooth. 

The teeth are treated as cylinders of lengths equal to the face 
width and radii equal to the radii of curvature at the contact point. 
For involute teeth the radii are the distances from the contact point 
to the tangent point of the respective base circles. For noninvolute 
profiles instantaneous base circles must be used. Fig\ire B3-6 depicts 
the previously discussed geometry considerations as applied to the ISG 
drive. Distances h^ and h^ are along the line of action from the con- 
tact point to the center of the teeth. 

Using Hertz's formulation for contact between cylinders 

- 8 Q r(l - y^)/irE . . . (B3-28) 

For the external-interned, gear combination 


1 1 



1 

r 


... (B3-29) 
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and 

Q = load per tooth face width 

= radius of curvatvire of external gear 
r^ = radius of curvature of internal gear 


also, 

max irb 

Considering the contacting geeu: teeth as slightly curved semi- 
infinite planes the Hertzian deformation 


E I 


Q 4(1 - u*) 
E ir 


[la 




2r(l - u )P 


1 - 


max 


. .. (B3-30) 


where 

6 and 6 are the deflection due to the Hertzian 
E I 

deformation of the external and internal tooth 
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Tooth Model for Hertzian 
Contact Deformation 



135 


APPENDIX C 

COMPUTER PROGRAM PACKAGE 


This section contains the computer listing of all three modules, 
typical output data and instructions for entering the data. 

C-1 Listing eind Sample Rm of the Static Analysis Program 
"Internal Static" 

C-2 Listing emd Sample Run of the Dyncunic Analysis Program 
"Interned Dyneunic" 

C>3 Listing emd Semple Run of the Stress Analysis "Internal 
Stress" 

C-4 Entering of Input Data 
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C-1 LISTING AND SAMPLE RUN OF THE STATIC ANALYSIS PROGRAM 
"INTERNAL STATIC" 
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IF (V2(L-1 ).GT.V1EP) GO TO 62 

IF (V2(L).LE.V1EP.AND.V2(L-1).LT.V1EP) GO TO 61 
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Y2EE=Y2(L2) 

I2SS=1 

I2EE=L2 



CALL DEFL(1,Y1EE,Y1SS, MEE. nSS) 
NNN=6 

CALL DEFL(2,Y2SS,Y2EE, I2SS, I2EE) 
NNN=9 
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T0EFL2(3,N)=RESIDL(2,N) 
8993 CONTINUE 

DO 901 J=1, IR1S22 
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or TOOTH ENGAGEMENT. IN EACH OF THESE TABLES POSITION 1 CORRESPONDS TO THE STARTING POINT OF CONTACT WHILE 
POSITION 50 CORRESPONDS TO THE END POINT OF CONTACT. 
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C-2 LISTING AND SAMPLE RUN OF THE DYNAMIC ANALYSIS PROGRAM 
"INTERNAL DYNAMIC" 
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SUBROUTINE V I BS( A, B, X, E IGV) 

IMPLICIT REAL*8(A-H,0-Z) 

ABS(X) = OABS(X) 

SQRT(X) = DSQRT(X) 

D I MENS I ON A( 4, 4 ) , B( 4, 4 ) , X( 4, *1 ) , E I GV( 4 ) , 0( 4 ) 
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IF(NP.EQ.4) NRK 

RETURN 

DT = DT/2. 


SUBROUTINE MORERK(X,DX,DT) 
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IF(T.GT.(PT-OT/2. )) GO TO 2 
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16 CONTINUE 

XTKMIN = .95*XTKMIN 


XTKMAX = 1.05*XTKMAX 
XTFMAX = 1.25«XTFMAX 
RK = XTKHAX - XTKHIN 
DELTAK = RK/4. 

RF = XTFMAX - XTFHIN 
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LINEK( I) = . 
155 LINEF( I )=JP 


WRITE{6, 1i|1 ) LIN£K,LINEF 
DO 165 I = l,l|l 
LINEK( I ) = JBLANK 
165 LINEF( I )=JBLANK 
RETURN 
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14.5 DEGREE PRESSURE ANGLE 
DIAMETRAL PITCH IS 8.000 
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THE SYSTEM USED IN THE DYNAMIC ANALYSIS IS PICTURED BELOV/. 
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THE SYMBOLS USED IN THE ABOVE FIGURE HAVE THE FOLLOWING DEFINITIONS: 



JC1 IS THE MASS MOMENT OF INERTIA OF THE DRIVING GEAR (GEAR 1) 

JGa IS THE MASS MOMENT OF INERTIA OF THE DRIVEN GEAR (GEAR 2) 

JL IS THE MASS MOMENT OF INERTIA OF THE ’LOAD 

KDS IS THE TORSIONAL SPRING STIFFNESS OF THE DRIVING SHAFT 

KG IS THE LINEAR SPRING STIFFNESS OF THE GEAR PAIR 
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THE NUMERICAL INTEGRATION WAS CARRIED OUT FOR A LENGTH OF TIME EQUIVALENT TO 140 CYCLES 
THE TIME REQUIRED FOR THE START-UP TRANSIENT TO DECAY (THIS TIME IS ASSUMED TO BE EQUALNG 
TO 5 TIMES THE LONGEST SYSTEM NATURAL PERIOD) TO THE TIME REQUIRED FOR ONE ADDITIONAL 
TOOTH PASSAGE CYCLE. THE DATA TABULATED IN TABLES 7 AND 8 BELOW COMES FROM THIS LAST 
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C-4 ENTERING OF INPUT DATA 

The three modules of the ISG drive computer package require input 
data vdiich are obtainable from the namelist arrays or the two data 
files. Numerical data may be entered without format statements, and 
fields are generated as required. The variables required for the name- 
list arrays along with their respective headings are: 

/HEDING/ 

TITLE 1 any suitable title or information 
TITLE 2 can be entered on three lines from 
TITLE 3 Title 1 through Title 3 

TAPE - Alphanumeric code to indicate whether certain 
information is to be filed on tape. Tape 8 is 
used for the dyneutiic analysis (Module 2) and 
Tape 9 is for the stress analysis (Module 3) . 

* YES * - write on tape 

*N0* - do not write on tape 

/CONTRL/ 

INPUT - alphanumeric code used to designate type of 
input data 

*ENGL' -English (Ibf, in., sec.) 

*SI* - metric (newtons, mrtf, sec.) 

OUTPUT - alphanumeric code used to designate output; 

codes used are same as for input 
MODE - alphaniameric code used to designate whether 
or not profile modifications are input 
*NO* - no modifications 

'YES* - modifications listed under /PRFDEF/ 
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/PHYPAR/ (two data points required per variable) 

E - Young * s modulus 

PR - Poisson’s ratio 

GAMA - specific weight 

JG - polar moment of inertia; optional, program 

will generate if no value entered 

/GENPAR/ 

DP - diametral pitch (English input only) 

M - gear module (metric module only) 

DELTP - backlash 

TIN - input torque 

RPMIN - input RPM 

ZETAS - damping coefficient of shaft 

ZETAG - damping coefficient between gear teeth 

PHID - pressxire angle (degrees) 

* JD " mass moment of inertia of driver 

* JL - mass moment of inertia of load 

* KDS - torsional spring stiffness of driving shaft 

* KLS - torsional spring stiffness of load shaft 

* LDS - length of drive shaft 

* LLS - length of load shaft 

/GEOPAR/ (two data points required per variable) 

TG - number of gear teeth 

AD - addendum 

WD - working depth 

GRRF - fillet radius of basic rack 

* RI - hub radius 
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FW - face width 

UCUT - undercut 

RT - rim thiclcness 

RADEL - radial deflection 

C0R12 - modifier for tooth center angle 

OOR34 - modifier for tooth center euigle 

/PARAME/ 

NLIM - cuigular sweep parameter - Gear 1 

MLIM - angular sweep parameter - Gear 2 

DELT - increment 

JJJJ - adjustaUsle do loop parameter - Gear 1 

T i T i T i T . - adjustcdsle do loop parameter - Gear 2 

DPSLIl - £uigular correction due to radial deflection - 
Gear 1 

DPSLI2 - 6uigxilar correction due to radial deflection - 
Gear 2 

DPELl - modifier for tooth center angle - Gear 1 

DPEL2 - modifier for tooth center angle - Gear 2 

* option^d./ program will generate if no value entered 

In addition to ev^d.uating purely involute gear teeth, the gear 
tooth profile can be modified to simulate tip relief or undercutting. 
Also, sinusoidal errors can be introduced, as well as pits, to simu- 
late involute errors due to manufacturing and surface damage, respec- 
tively. These modifications are introduced in the /PRFDEF/ namelist. 
If MODF » NO, /PRFDEF/ is not included in the data card set. 
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/PRFDEF/ (two data points required per variable) 

PATH - parabolic tip modification 

STTM - straight line tip modification 

RATM - roll angle of tip modification 

PABM - parabolic bottom modification 

STEM - straight line bottom modification 

RABM - roll angle of bottom modification 

PER - amplitude of sinusoidal error 

PAP - phase angle of sinusoidal error 

CYC - number of cycles of sinusoidal errors 

IPrr - profile coordinate points over which pit 

occurs 

DEEP - depth of pit 

Use of the namelist arrays offers a simple, unformatted means of 
inputting data and is very convenient for looping more than one data 
set. After the initial data set, siibsequent data sets need just two 
input revisions. If a later namelist array contains no revisions, 
only a card with the eurray heading and ending need be submitted. Un- 
listed variables default to the previous values. Examples of input 
data card sets illustrate the following namelist data card format. 

1. Column one is blank. 

2. is used to signify new naunelist array* 

3. is followed by the namelist name. 

4. A blank separates the namelist name and the first 
vauriaible name. Subsequent variaUdles are separated 
by commas. 
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There are two methods for defining the two element 
variables. The elements are defined in the order they 
are to be entered in the variable and separated by 
commas, i.e. , TG=32 , 96 defines TG(1) = 32 and TG{2)=96. 

If both elements are equal, they may be entered by 
listing the number of identical values, the multipli- 
cation symbol, and then the value itself, i.e., 

AD=2 *0.125, defines AD(1) =0.125 and AD (2) =0.12 5. 

6. The last listed array value is followed by a blank and 
then the symbol from column 2 is repeated. The word 
END immediately follows the symbol and signifies the 
end of that array. 

Because of the modular approach it is necessary to store certain 
information from one module for use in another module. Tape 8 stores 
the pertinent data from the static auialysis which is needed for the 
dynamic analysis. Tape 9 stores the pertinent data from the dynamic 
analysis for use in the stress analysis. This tape contains some of the 
previously transferred static anadysis data on Tape 8. The stress analy- 
sis then has logic to initiate a static or dynamic stress analysis. 

The modules have the capability for accepting either SI or English 
gear input data and have options to print the results in either SI or 
English units. Input and output do not necessarily have to be of the 
same regime, i.e., SI output can be obtained from English input and 
vice-versa. Data submitted under the 'ENGL* code should be in pounds- 
force, inches and seconds. The data submitted under the *SI* code 
should be in newtons, millimeters and seconds. The only exception to 
this is the density value under the *SI* code should be in kg/m . 
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GLOSSARY OF TERMS 


Text 

Computer Prograun 

a 

AD 

A 


b 


B 

BGM 

c 

CR 


C 


Cb 


CYC 

d 

DED 

DFl, DF2 

DPI, DF2 

D 


deep 

DEEP 

E 

E 

F 

FW 

F 


r 


F 


t 


G 


h 


HSF 

HSF 

K 

HT 

t 


I 


ISG 


J 


[J] 

[MM] 


KG 

KP 


Description 

addendtm 

area 

ring gear width 

maximum thickness of tooth 

backlash 

loaded contact ratio 

center distcuice of gears 

damping coefficient 

number of sinusoidal error cycles 

dedendum 

dynamic load factors 

ball diameter 
depth of pit 

Young's modxilus 

face width 
geometry factor 

radial load of bearings 

tcuigential load of gears 

modulus of torsion 

ring gear thickness 

hub/ring torsional stiffness factor 

whole depth 

moment of inertia 

internal gear drive 

polar mass moment of inertia 

inertia matrix 

gear mesh stiffness 

gear pair stiffness 
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Text 

Computer Program 

Description 

[K] 

[SM] 

stiffness matrix 

L 


length of roller 

L.A. 


line of action 


MG 

gear ratio 


CR 

contact ratio 

M 


bending moment 
module 

N 


normal load 


TG 

number of teeth 


OMG 

constant angrilar velocity 

P 

CP 

circular pitch 

Pb 

BP 

base pitch 

P 


applied load at contacting point 

PABM 

magnitude of parabolic modification 
at bottom 


PAP 

angle from start of sinusoidal error 
to start of involute 


PATH 

magnitude of paraO^olic modification 
at tip 


PE 

profile error 


PER 

maximxim profile error 


PH 

Hertzian pressure 


PM 

profile modification 


PPD* 

instantaneous pitch radius 


PSITP 

static cuigular position 

r 

RPC 

radius at pitch circle 

^A 

RAC 

radius at addendum circle 

r. 

RBC 

base radius circle 

b 

RABI 

RA at the bottom of involute 


RABM 

length of root modification in 
degrees of roll 


RABOT 

RA at bottom of involute 


RAM 

roll angle at end of modification 
at tip 


RAN 

RA at end of modification at bottom 


RAPP 

RA at pitch point 


RAT 

RA at tip 


RAT IP 

RA from end of modification to RA 


at the pitch point 

length of tip modification in degrees 
of roll 


RATM 
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Text 

Computer Proqram 

Description 

RCCP 

RCC 

instantaneous radius of curvature 


RF 

fillet radius 


RLM 

limit radius 


RT 

edge radius of generating tool 

Q 

Q(k). 

static gear pair load 

2 t 

QT 

total mesh static load 

QD 

QD(k)^ 

dynamic gear pair load 

QDt 

QDTd) 

total mesh dynamic load 

®max 


maximum radial bearing load 


STEM 

magnitude of straight line modifi- 
cation at bottom 


STTM 

magnitude of straight line modifi- 
cation at tip 

SV 


sliding velocity 

T 


torque 

T 

IN 

TD 

input torque 

T 

LA 


theoretical line of action 

TR, TR' 

TDIN 

theoretical and instantaneous 
transmission ratio 

’’'out 

TOUT 

output torque 

Ui, 02 


interval of contact 

u 


cdDscissa in global coordinate 
system 

V 


shear 

velocity 

ordinate in globad coordinate 
system 

w 


abscissa in rotating coordinate 
system of gears 

X 


abscissa in local tooth coordinate 
system 



Conipvtter Program 


Description 


ordinate in locad tooth coordinate 
system 

ordinate in rotating coordinate sys 
tern of gears 

number of balls or rollers-bearings 


GREEK SYMBOLS 


TDEFL 

RA 

PSP 

PSPD 

PSPDD 

PR 

ZETAG 

ZETAS 

PI 


contact angle - bearings 

angle between point of contact emd 
center of tooth 

angle between pitch point and center 
of tooth 

maximum fillet stress angle 

deflection 

roll angle 

involute polar angle 

dyncunic displacement 

dyncutiic velocity 

dynamic acceleration 

Poisson's ratio 

critical damping ratio - gear 

critical damping ratio - shafts 

3.141592654 


stress 

torsional stress 
pressure angle at any point 
normal pressure angle 
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Text 

1 

2 

D 

G 

i 

k 

s 



external gear 
internal gear 
driving element 
gear 

mesh arc position 
tooth number 
shafting 
instanteuieoxis 
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